NASA Contractor Report 182256 


1 6. 

/ / A 9 ^ ^ 

/ 


Experiment Plans to Study Preignition 
Processes of a Pool Fire 
in Low Gravity 

»69- 19442 

Dnclas 
0195460 


Lrf7f.T; 1fc ^ 25ej iXPIBlHEBl ■ iIAS6 10 SltJDl 

F !? C M £££ Cf A £CC1 FJEE 1B W* 

e.6. lfc €£J£ - 1S66 fiE2j Be port 
(Cdiitoirid Oiiv.) 1i6 p CSC1> 22A 


G3/29 


David N. Schiller 

University of California— Irvine 
Irvine, California 


March 1989 


Prepared for 

Lewis Research Center 

Under Grant NAG3-627 


NASA 

National Aeronautics and 
Space Administration 


Contents 


Table of Symbols 

Acknowledgements iv 

Chapter 1 Introduction 1 

Chapter 2 Summary of the Science 3 

2.1 Background 3 

2.2 Literature Review 8 

2.3 Justification for Space-Based Experiment 15 

2.4 Objectives 16 

2.5 Use of NASA’s Microgravity Facilities 17 

Chapter 3 Identification of Experiments 18 

3.1 Test Matrix 18 

3.2 Measurement Requirements 24 

3.3 Experimental Conditions 36 

3.4 Preliminary Hardware Concepts 40 

Chapter 4 Experimental Timeline 41 

4.1 1-g Tests 41 

4.2 Drop Tower Tests 42 

4.3 Space-Based Tests 43 

Chapter 5 Planned Data Analysis 44 

Chapter 6 The Numerical Model 45 

6.1 Description of Numerical Model 45 

6.2 Governing Equations 49 

6.3 Source Term Approximations 59 



Chapter 7 Preliminary Results 66 

Chapter 8 Conclusions 81 

Appendix A Description of Preliminary Hardware Concepts .... 83 

Appendix B Derivation of Energy Equation 101 

Appendix C Derivation of Source Term For Species Equation With 
One-Step Chemical Reaction 117 

References 120 


ii 


Table of Symbols 


Gr Grashof number. 

Ma Marangoni number. 

Re Reynolds number. 

Pr Prandtl number. 

Eu Euler number. 

Ec Eckert number. 
fi Volume expansion coefficient. 

g Gravitational acceleration (or effective acceleration with 
respect to an accelerating frame of reference). 

g Gravitational acceleration at sea level (9.81 m/s 2 ). 

L Characteristic length (typically taken to be the radius of 
the test cell). 

Tf Maximum fuel temperature at the liquid surface. 

T 0 Bulk temperature of the fluid. 

Th Heater temperature. 
p Density of the fluid. 
v Thermal diffusivity of the fluid, 
a Kinematic viscosity of the fluid. 

a Surface-tension coefficient. 

U Reference velocity. 

6 Local angle of the liquid surface. 

8 Liquid surface deformation with respect to the height of 
the liquid in contact with the side walls. 

A T Temperature difference — T a . 

Note : Many other symbols are used and defined in Chapter 6. 
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Introduction 


The purpose of this document is to specify requirements to guide experimental 
studies of the fundamental behavior of a liquid pool fire prior to ignition. These 
specifications are provided for experiments in three separate settings: (a) a normal 
gravity laboratory; (b) the NASA-LeRC drop towers, and (c) a space-based laboratory 
(e.g., Shuttle, Space Station). A rationale is developed for both minimum and desired 
requirement levels. In general the requirements are most stringent in the normal 
gravity experiments. 

In Chapter 2, the scientific problem is described with the intention of justify- 
ing the need for low gravity experimentation. Chapter 3 identifies the test matrix 
and measurement requirements for each experimental setting. Chapter 4 provides a 
typical operational sequence for the experiments. The planned data analysis is de- 
scribed in Chapter 5. Preliminary hardware concepts are described in Appendix A. 
The numerical model is described in Chapter 6, which also includes the governing 
equations and source term approximations used in the code. Preliminary numerical 
and experimental results are then given in Chapter 7. 
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All experiments in this document are intended to collect scientific data. The re- 
quirements specified herein will be reviewed and updated periodically (e.g.. annually) 
as appropriate. 
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Chapter 2 

Summary of the Science 


2.1 Background 

Fire safety problems of ignition and flame spread above a pool of liquid fuel 
arise in many accident situations where liquid fuel spills in the vicinity of hot ignition 
sources. Another prototypical hazard situation might be the rupture of a fuel tank 
when hot engine parts or exhaust gases may appear in close proximity of liquid fuel. 
Heat and mass transport which control ignition delay and flame spread rates are 
complicated by multiple energy and mass transport processes, phase change, and 
chemical reaction. It is believed that reduced gravity will remove or isolate some of 
these complications. However, the scientific literature is devoid of experimental study 
of these characteristics under conditions of reduced gravity. 

The literature of liquid pool burning theory has for the most part neglected 
the gas and liquid phase coupling, and concerned itself solely with liquid motion. 
However, the most interesting scientific situation (and one of very practical impor- 
tance) involves not-so- volatile fuels whereby the coupling between the liquid and gas 
phases is of critical importance. That is, heat from the gas phase is required to 
vaporize some portion of the liquid fuel in order to provide a combustible mixture 
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in the gas phase. Convective heat and mass transport between the liquid and gas 
phases is critical and often dominant during the preignition period. Other modes 
of heat transfer, such as thermal radiation from the ignition source, may also play 
a very important role in the liquid heating. Temperature gradients in the gas and 
liquid phases, as well as concentration gradients in the gas mixture, create natural 
convection currents. Temperature gradients in the liquid surface initiate “surface- 
tension-driven convection” which can be a critical factor in the ignition delay and 
flame spread phenomenon. On account of nonuniform heating from above due to 
the ignition source or the flame, liquid motion will be driven by both surface-tension 
gradients and by buoyancy. These driving forces generally support in concert surface 
fluid motion away from the heat source. Therefore this motion tends to delay ignition 
as heat from the source is converted rather than concentrated. On the other hand, 
the motion supports flame spread as the convection assists the preheating process 
ahead of the flame. The gravity effect on ignition and flame spread is, in general, 
opposite to that of surface tension. As gravity decreases, the fluid mixing due to 
natural convection in both phases is retarded and the liquid at the surface is heated 
and vaporized more rapidly. Therefore the reduction in gravity should lead to shorter 
ignition delay and slower flame spread. However, the final result depends on the rela- 
tive strengths of buoyancy versus surface tension. Surface-tension effects will tend to 
dominate over buoyancy at reduced gravity and smaller characteristic lengths. It is 
obvious, therefore, that the magnitude of the gravitational acceleration is one of the 
most important parameters that can be studied. 

In order to understand the physical phenomena involved in ignition arid flame 
spread, W. Sirignano et al., under NASA/OSSA sponsorship, developed a compu- 
tational code which has been used to study the preignition motion of a deep liquid 
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fuel pool in an enclosed container heated unevenly from above. The key findings of 
the research are contained in References [1] and [2]. If the effect of radiation is not 
considered, then as the Grashof number decreases (or as the gravity force diminishes) 
the number of circulatory cells and the maximum value of the velocity in both the 
liquid and the gas phase decrease markedly. The convective heat transfer dominates 
conduction and the buoyancy-driving of the flow is greater than the surface-tension 
effect at high Grashof numbers (higher than 10 6 ). At lower Grashof numbers (10 3 
and less) conduction is dominant. Surface tension dominates buoyancy in inducing 
liquid motion at lower Grashof numbers, but this motion is itself not as important 
as gas phase conduction. These conclusions apply to pools of depths comparable to 
their widths. It is expected that surface tension will be more important in shallow 
pools. Radiation and thermal inertia of the walls were shown to produce significant 
quantitative modifications. An important qualitative difference is that as the con- 
tainer walls are heated by radiation from the primary heat source, strong natural 
convection currents are initiated near these hot surfaces. The circulatory cells near- 
est the primary source need not be the strongest vortical structures when buoyancy 
is present. 

The geometry of the problem is shown on Figure 2.1. The initially quiescent 
fluid in the container consists of two phases: liquid and gas. A constant-temperature 1 
heater surrounded by a water-cooled jacket provides nonuniform heating from above. 
Once heating begins, the response of the system in terms of velocity and temperature 
is simulated at 400 node points (or 1600 node points using the CRAY computer) in 

both the gas and liquid phases. The effects of gravity are accounted for by varying 
'The heater actually need not be exactly constant temperature as long as the radial temperature 
profile of the heater is known. 
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Figure 2.1: Geometry of the problem (taken from Reference [1]). 


the Grashof number. Because the prediction of the liquid-vapor interfacia] shape is 
complex, the model presently considers only a flat liquid-vapor interface. Thus, as 
discussed later, provisions must be made in designing reduced gravity experiments to 
maintain a flat liquid- vapor interface. 


The general philosophy of the computational effort has been to simulate the 
more tractable aspects of the problem initially, and to increase the complexity of the 
computer program in small steps to approach the real case. First, heat transport 
in a one-phase fluid within an enclosure was studied using the Los Alamos SOLA 
code [2]. This code was modified to account for natural convection. Then, tran- 
sient natural convection in a two-layer gas-and-liquid system without surface tension 





was studied. Next, surface-tension effects were added. In the fourth iteration, ra- 
diation was added. Inclusion of radiation effects caused the necessity of revision of 
the thermal boundary condition at the walls of the container. Thus, the effects of 
thermal inertia and one-dimensional conduction along the walls were added. The 
following steps in the research were to include the effects of variable properties in the 
gas and liquid phases, evaporation and mass transfer, and two-dimensional thermal 
conduction along the walls. This task, however, required a complete change of the 
numerical approach and the computer code. The incompressible SOLA code could 
not be employed to solve the variable property fluid phenomena. A new computer 
code therefore was developed based on the SIMPLE algorithm which originated at 
London Imperial College [3]. This computer program has a capability to solve tran- 
sient, multicomponent flows at low Mach numbers including the effects of variable 
thermophysical properties (density, viscosity, etc.). Variable density is important in 
the gas phase because during preignition heating, high temperature gradients are ex- 
pected near the ignition source (i.e., the heater). In such a situation, the use of the 
constant density approach along with the Boussinesq approximation may produce 
very serious quantitative errors in the analysis of the phenomenon. The incorpora- 
tion of variable thermophvsical properties in the liquid phase is also very important 
since the viscosity and volume expansion coefficient of many liquids vary considerably 
with temperature. The numerical method used in the current computational model 
is described in Chapter 6. In the following section a literature review of ignition and 
flame spread at normal and reduced gravity conditions is presented. 
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2.2 Literature Review 


Comprehensive reviews of the state-of-the-art on flame spreading across liquid 
fuels have been made by Glassman and Dryer [4] and Kanury [5]. Systematic experi- 
mental studies of the ignition phenomenon were undertaken by Princeton researchers 
in the seventies [6,7,8]- It was found that the factors affecting the ignition and the 
flame spread depend strongly on whether the bulk fuel temperature is above or below 
its flash point. Conceptually, the flash point is the liquid temperature which creates a 
lean flammability limit mixture of fuel vapor and air over the liquid surface. The flash 
point can be calculated for pure fuels if the lean flammability limit and the saturated 
vapor pressure are known as a function of temperature. However, the experiments 
with open pools generally give the flash point temperature several degrees 

higher than the theoretical value. The experimental values of Tfj^^ are strongly 
affected by the height of the ignition source above the liquid surface. 

When the liquid fuel temperature is below r the flash point, the ignition delay 
depends upon many factors, such as the initial liquid temperature, liquid viscosity, 
depth of the pool, and surface tension of the liquid. The heat flux from the igni- 
tion source to the liquid is usually maximum just beneath the source and decreases 
monatonically with the distance from the igniter. The nonuniform heating causes a 
temperature gradient along the fuel surface. Since the surface tension of liquid varies 
inversely with temperature, the gradient of surface tension will induce a surface flow 
from the warmer temperature to the cooler temperatures. Experimentally observed 
surface velocities have been reported of the order of several centimeters per second. 
The hot liquid fuel below the igniter is continuously carried away by fluid motion 
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and replaced by the cold liquid from the volume of the pool. Therefore, the surface- 
tension-driven convection results in a considerable retardation of the ignition process. 
Suppression of the motion by artificially reducing surface tension or increasing liquid 
viscosity enhances the ignitability markedly [4,6]. 


In contrast with the ignition phenomenon, flame spread across the pool is as- 
sisted by the surface-tension-driven convection which can be a dominant factor con- 
trolling the flame spreading rate [4,8]. In this case, the flame above the pool plays the 
role of the ignition source. Hot liquid moves ahead of the flame front thus promoting 
the ignition of the fuel vapor. Theoretical studies of steady-state flame propagation 
controlled surface-tension flows have been made by Sirignano and Glassman [8] and 
Torrance [9,10]. Torrance considered also the influence of natural convective motion 
in the liquid induced by nonuniform heating from above. The buoyancy effects inter- 
act with the surface-tension flows. Furthermore, in flows where the Reynolds number 
is greater than approximately 10 2 or 10 3 , inertial forces balance the buoyancy force. 
The relative importance of buoyancy, surface-tension, and inertial effects can be an- 
alyzed using the nondimensional Grashof (Gr), Marangoni (Ma), and Reynolds (Re) 
numbers: 


Gr 

Ma 

Re 


(WAT) - Tp) 

- (h - (?) 

pua \aT ) \ax ) 3 

UL 


Here L is the characteristic dimension of the pool; p, y, a are the density, kinematic 
viscosity and thermal diffusivity of the liquid, respectively; T 0 is the bulk liquid 
temperature; Tj is the maximum fuel temperature at the liquid surface (may be 
replaced by the heater temperature. T),, when the effect of radiation is included); a 
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is the surface-tension coefficient,; ft is the volume expansion coefficient: and V is a 
reference velocity. Note that the Grashof and Reynolds numbers may similarly be 
defined for the gas phase. In this case gas phase properties are used, Tj is replaced 
by T h , and T 0 denotes the bulk gas phase temperature. The nondimensional vector 
form of the Navier-Stokes equations using the Boussinesq approximation for scaling 
purposes may be written as 

— + (u -V)u = — (?-)t - EuVP + -}-V 2 u 
dt v ’ Re 2 \g) Re 

Thus when inertial forces balance the buoyancy force one obtains 

Gr ~ Re 2 

The order of magnitude of the liquid surface velocity caused by the buoyancy force 
may therefore be estimated as 

U Te f ~ t'VGr/L - \fpgL\T 

Alternatively, a reference velocity may be derived by balancing the tangential shear 
stresses with the surface-tension gradients at the free surface. For /?e > 1 a boundary- 
layer flow will occur and the resulting reference velocity, U„, may be estimated as [II] 

rr v(Ma/Pr) 2 / 3 

V ° j 

It should be noted, however, that scaling of the phenomenon is a difficult task because 
the values of 7/ and ( dT/dx) s are not generally known a priori. This is especially true 
for the ignition problem which is transient by its nature and where the Grashof and 
Marangoni number vary with time. For this reason, the above reference velocities 
cannot be used in the computational model. However if T/ is replaced by Th in 
the definition of the Grashof number, then (7 re f may be estimated a priori by its 
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above definition. A very rough estimation of Gr for liquid R113 in 10 -4 g with the 
arbitrarily chosen values of L = 10 cm and Th~T 0 = 500°C yields: Gr = 1 x 10 5 and 
U Te { ~ 2 mm/s. Torrance [9] found that in a liquid-only simulation at Gr ~ 10 6 , the 
buoyancy influence on the liquid flow patterns in the subsurface layer is comparable 
to the surface-tension effect. 

Recently, Furuta et al. [12] performed a finite-difference study of a steady lami- 
nar flame spread over a shallow fuel pool. The momentum, energy, and mass conser- 
vation equations were solved simultaneously both in the gas and liquid phases. The 
combustion process was modelled using the flame-sheet approximation. The results 
of the calculations showed that buoyancy and volume expansion effects influence the 
gas flow in the region ahead of the flame front. The surface-tension effect on the 
interface induces recirculation motion in both liquid and gas phases. Such a coupled 
circulation structure at the gas-liquid interface plays an important role in the flame 
spread process [4]. 

It should be noted that in contrast to steady-state flame spread, the transient 
ignition process has not been treated theoretically in the literature. The preigni- 
tion hydrodynamics and heat transfer phenomena have been numerically studied by 
Sirignano et al. in work sponsored by NASA [1]. Detatiled descriptions of this study 
are given in [1] while the most important results are indicated below. 

The calculations were made for the cylindrical enclosure (see Figure 2.1) con- 
taining a liquid-air two-layer system. Approximations made in the numerical model 
include the following: 

1. The Boussinesq approximation is made with constant thermophysical properties 

in both the gas and liquid phase. 
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2. The 1 mm mild steel walls bounding the enclosure are assumed to be so thin that 
the temperature drop across the thickness is negligible (the “fin” approxima- 
tion). The top and side walls are adiabatically isolated from the environment; 
however, the wall temperature may change spatially and with time as a result of 
the convective heat exchange with the fluids, the radiative transfer between the 
different surfaces, and the thermal conduction along the walls (one-dimensional 
transient heat conduction). 

3. All of the solid surfaces as well as the gas/liquid interface are assumed opaque 
and black. Gas radiation is neglected. 

4. Liquid vaporization, diffusion, chemical reaction, and radiative absorption in 
the gas phase are neglected. 

5. The liquid surface is assumed to be flat. 

The typical streamline and temperature patters are illustrated in Figures 2.2 and 2.3 
for a liquid decane-air system at normal gravity conditions. Here, the dimensions of 
the cylindrical enclosure are: R = 10 cm, H = 20 cm and height of the liquid is 10 cm. 
The enclosure is heated by a small circular source (r = 2 cm) at the center of the top 
cover. The hot spot has no thermal contact with the remaining part of the top cover. 
The initial temperature of the system is 300 K, while the temperature of the heat 
source is held constant at 1000 K. The bottom wall of the enclosure is maintained at 
the constant temperature 300 K. 

In the liquid phase, the major changes of velocity and temperature occur in 
a thin region very close to the surface. Below this boundary layer, the liquid mo- 
tion exhibits the multi-vortex structure which is caused by the interaction between 
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Rodius,cm 


Figure 2.2: Typical streamline patterns for a liquid decane-air system at normal 
gravity conditions (taken from Reference [1]). Sign convention for nondimensional 
stream function shown on figure is opposite to equation (7.1). 
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Figure 2.3: Typical temperature contours for a liquid decane-air system at normal 
gravity conditions (taken from Reference [1]). 



the surface-tension-driven and natural convection phenomena. Contrastingly, only 
single vortices are observed in the gas and liquid phases at low gravity conditions 
(g = 10 -4 m/s 2 ). Figure 2.3 shows that the hot gas zone is concentrated in close 
proximity to the ignition source. Below this zone the gas remains relatively cold. A 
slight temperature rise is also observed near the liquid surface which is heated by 
radiation. Therefore, the ignition of fuel vapor/air mixture will probably be initiated 
near the hot source rather than close to the liquid fuel surface. 

The calculations demonstrated the great importance of including effects of ther- 
mal radiation and thermal inertia of the walls bounding the enclosure. Nonuniform 
radiative heating of a liquid surface appears to be the principal factor in the initiation 
of a surface-tension-driven convection along the gas-liquid interface. The thermal in- 
ertia of the walls controls the temperature growth in the gas phase. The gas and 
liquid motions are the net results of both surface-tension and buoyancy forces. 


2.3 Justification for Space-Based Experiment 

As mentioned in the previous section, the scientific literature is devoid of exper- 
imental study of ignition and flame spread characteristics under conditions of reduced 
gravity. It is not clear whether ignition susceptibility and flame speed will increase or 
decrease in reduced gravity. The only way to change the scale of experiments (particu- 
larly to decrease the Grashof number) without affecting other important parameters 
is to change the gravity level. Currently planned ground-based low gravity (drop 
tower) experiments, though necessarily transient in nature, should provide an excel- 
lent basis for assessing relative effects of various gravity levels and determining the 
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feasibility of proceeding to a space-based experiment. However, modelling suggests 
that reduced gravity experiments will require several minutes to reach a quasi-steady 
state. A quasi-steady state condition is characterized by stabilized thermal and dy- 
namic boundary layers [1]. For a 10 cm height of air at 300 K, the time scale for 
thermal diffusion may be estimated as 

1 2 

T ~ — = 454 seconds, 
a 

The significance of this is that it would take on the order of 7-8 minutes for heat to 
conduct through a 10 cm height of air. Therefore the full effect of heat conduction 
cannot be realized in a 5 second drop tower experiment. 


Quasi-steady state is important in both the ignition and flame spread problems. 
In the current preignition problem the initial liquid temperature is below its vapor- 
ization temperature and well below its flash point temperature. Hence a relatively 
long period of time (^> 5 seconds) is required to study the asymptotic behavior of 
the liquid fuel as it vaporizes. The quasi-steady state behavior of the experiment will 
give insight to the associated flame spread problem, where the primary interest is the 
asymptotic spreading rate of a flame once it has been established. In order to study 
the quasi-steady state behavior of the preignition phenomenon under the required 
range of 0gL&T and the Grashof number, a space-based facility as offered by the 
Shuttle is required. 


2.4 Objectives 

The principal objective of the reduced gravity experimental effort is to ver- 
ify experimentally the predicted behavior of an enclosed liquid fuel pool during the 
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preignition phase and to determine via measurement the role of gravity on such be- 
havior. Interesting questions that should be resolved by the research concern: (?) the 
determination of domains where conduction, convection, or radiation are dominant. 
( ?'?') the time required to reach an ignition, ( Hi ) the importance of heat transport 
through the gas versus heat transport through the liquid, and (iv) the domains where 
surface-tension-driven flow dominates buoyancy-driven flow. Space-based experimen- 
tation will be supported by the aforementioned computational research as well as 
ground-based experimentation. NASA-LeRC will conduct 1-g tests and drop tower 
tests to support the research program. By varying the majority of experimental pa- 
rameters in these ground-based test s, the cost of space-based experimentation will be 
minimized. The science requirements for both the space-based and the ground-based 
experimentation will be given in this document. 


2.5 Use of NASA’s Microgravity Facilities 

The NASA-LeRC 5.18 and 2.2 second drop tower facilities will be used to obtain 
low-gravity transient data to support the computational and space flight experimental 
efforts. Freon 113 (liquid) and air will be used in plexiglass containers of planar and 
cylindrical geometries. Data from the drop tower tests will be useful in verifying the 
design and feasibility of a space-based experiment. For example, the hydrodynamics 
of maintaining a flat gas-liquid interface in reduced gravity may be evaluated in drop 
tower tests. The effectiveness of flow visualization techniques may also be assessed 
in these short duration microgravity tests. However, the dynamics of the transient 
reduction of gravity in the drop tower tests may prevent valid scientific data collection 
in the available reduced gravity time. 
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Chapter 3 


Identification of Experiments 


3.1 Test Matrix 

The following section summarizes the test matrix for each of the three experi- 
mental settings (1-g, drop tower, and space-based). 

The minimum goal of the entire experimental research is to study the problem 
of interest over at least a three order of magnitude range in the inertia of the flow. 
When Re > 10 2 or 10 3 , inertial forces balance the buoyancy force and the inertia of 
the flow may be estimated by 

U 2 xe{ ~ ,3gL(T h - T 0 ) = dgLAT 

where all variables are as before except 7\, which is the heater temperature. This 
range should cover Gr = 1 0 3 — 1 0 6 in order to study the relative effects of convection 
versus conduction and buoyancy versus surface tension on the flow. The transitions of 
dominating flow effects cited earlier in Section 2.1, namely Gr ~ 10 3 and Gr ~ 10 6 , 
require experimental validation. Contrary to the previously given definition of Gr 
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(which used the maximum fuel temperature of the liquid surface as the hotter tem- 
perature), the Grashof number cited above is determined by 

^ / 1gL 3 (T h - T 0 ) 

v i 

It would be more desirable to cover a four or five order of magnitude range in 0gLAT 
to provide more assurance that the key ranges of buoyancy and inertial forces as well 
as the transitions are studied in the experiments. The characteristic length scale for 
j3gLAT and the Grashof number is somewhat arbitrary. In previous simulations of 
a cylindrical test cell, the radius of the container was chosen. In future experiments 
with a rectangular test cell, the hydraulic diameter or the half-length of the container 
may be chosen for the length scale. Whatever length scale is chosen, it is the quantity 
/ 3gL\T (= > n combination with other important parameters (e.g., aspect ratio, 

heater size, etc.) which drive the results. Although the drop tower provides about a 
five order of magnitude reduction in gravity (and hence 0gLAT), the results of such 
tests are only transient. Through the use of the Shuttle, one can obtain a reduction 
of gravity of three or four orders of magnitude. Varying the heater temperature 
(AT — Th, — T 0 ) will provide an additional order of magnitude variation. Changing L 
will also linearly affect the inertia of the flow. The test matrix given in Table 1 was 
designed to (1) provide for the desired variation in flgLAT for the quasi-steady state 
experiments, and (2) investigate issues and effects of parametric variations relevant 
to the design of the space- based experiment. 


I 
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Table 1A: Test Matrix 


Variable 

1-9 

Drop Tower 

Space- Based 

Test Cell 

Cylindrical 

Cylindrical 

Cylindrical 

Dimensions 

(H. x dia.): 

(H. x dia.): 

(H. x dia.): 


• 1 cm x 1 cm 

• 1 cm x 1cm 

• TBD 


• 5 cm x 5 cm 

• 5 cm x 5 cm 

(10 cm x 10 cm) 


• 10 cm x 10 cm 

• 10 cm x 10 cm 



Rectangular 3 

Rectangular 3 (single phase only) 


Note: L is half-length 

(H x 2L x W): 

(H x 2L x W): 


of test cell (see 

• 20 cm x 15 cm 

• 10 cm x 7.5 cm 


Figure 3.2) 

x 30 cm 

x 15 cm 



• TBD 

• TBD 


Test Cell 

• Plexiglass, 1/8- 1/4” 

• Plexiglass, 1/8-1/4” 

• TBD 

Material 

thick 

thick 

(plexi glass, glass, 


• Thin glass 3 (TBD 


or sapphire) 


thickness) 



Liquid Phase 

• Silicon oil 

• Silicon oil 

R113 with glass or 

Fluid 

• R113 

• R113 

sapphire test cell 


• Water “ 



Gas Phase Fluid 

Air 

Air 

Air 

Heater Dimensions 

Cylindrical (dia.): 

Cylindrical (dia.): 

Cylindrical (dia.): 

(not including 

• 5 mm 

• 5 mm 

• 2.5 cm 

water-cooled 

• 1 cm 

• 1 cm 

• 5 cm 

jacket) 

• 2.5 cm 

• 2.5 cm 



• 5 cm 

• 5 cm 



Rectangular (length) 6 : 

Rectangular (length) 6 : 



• 2.5 cm 

• 2.5 cm 



• 5 cm 

• 5 cm 


Heater 

• 350 K 

• 350 K 

• 350 K 

Temperature 

• 575 K 

• 575 K 

• 800 K 


• 800 K 

• 800 K 


Heater 

• 1.0 

• 1.0 

• 1.0 

Emissivity 

• TBD 



Liquid Fill 

• 0.125 

• 0.125 

• 0.125 

Level (fraction 

• 0.25 

• 0.25 

• 0.5 

of container height) 

• 0.5 

• 0.5 


Use of Tracers 

• Yes 

• Yes 

• Yes 

for Flow 

• No 

• No 


Visualization 




Use of Thermo- 

• Yes 

• Yes 

• No 

couples for 

• No 

• No 


Measurement of 




Fluid Temperatures 





“Optional. 

6 The rectangular heater strip with water-cooled jacket spans the width of the rectangular container. 
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Table IB: Sequential Test Plan 


A. Single-Phase, 1-g Tests 


Test Celt 1 

Heater Size 

Heater Temp. 

Gas Tracers 

• 10 cm x 10 cm 

• 1 cm 

• 2.5 cm 

• 5 cm 

• 800 K 

• 575 K 

• 350 K 

Tobacco smoke 

• 5 cm x 5 cm 

• 1 cm 

• 2.5 cm 

• 800 K 

• 575 K 

• 350 K 

• 1 cm x 1 cm 

• 5 mm 

• 575 K 

• 350 K 


"Test cell is cylindrical, plexiglass unless otherwise specified; dimensions given are Height x Diameter; 
additional tests with rectangular cells and/or thin glass material are optional. 


B. Two-Phase, 1-g Tests 


Test Celt 1 

Heater Size 

Heater Temp. 

mrnmm 

Liquid Fill LeveT 

Gas/ Vapor Tracers 

• 10 cm x 10 cm 

• 1 cm 

• 2.5 cm 

• 5 cm 

• 800 K 

• 575 K 

• 350 K 


• 0.5 

• 0.25 

• 0.125 

• None 

• Tobacco smoke 

• 10 cm x 10 cm 

• 1 cm 

• 2.5 cm 

• 5 cm 

• 800 K 

• 575 K 

• 350 K 

R113 

• 0.5 

• 0.25 

• 0.125 

• None 

• Tobacco smoke 

• R113 vapor 

• 5 cm x 5 cm 

• 1 cm 

• 2.5 cm 

• 800 K 

• 575 K 

• 350 K 

Silicon Oil 

• 0.5 

• 0.125 

• None 

• Tobacco smoke 

• 5 cm x 5 cm 

• 1 cm 

• 2.5 cm 

• 800 K 

• 575 K 

• 350 K 

R113 

• 0.5 

• 0.125 

• None 

• Tobacco smoke 

• R113 vapor 

• 1 cm x 1 cm 

• 5 mm 

• 575 K 

• 350 K 

Silicon Oil 

• 0.5 

• None 

• Tobacco smoke 

• 1 cm x 1 cm 

• 5 mm 

• 575 K 

• 350 K | 

R113 

• 0.5 

• None 

• Tobacco smoke 


"Test cell is cylindrical, plexiglass unless otherwise specified; dimensions given are Height x Diameter; 
additional tests with rectangular cells and/or thin glass material are optional. 

Additional tests with water as the liquid are optional. 
c Liquid fill level is expressed as fraction of container height. 
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Table IB: Sequential Test Plan (continued) 


C. Single-Phase, Drop Tower Tests 


Test Cell a 

Heater Size 

Heater Temp. 

Gas Tracers 

• 10 cm x 10 cm 

9 1 cm 

• 2.5 cm 

• 5 cm 

• 800 K 

• 575 K 

• 350 K 

Tobacco smoke 

• 5 cm x 5 cm 

• 1 cm 

• 2.5 cm 

• 800 K 

• 575 K 

• 350 K 

• 1 cm x 1 cm 

• 5 mm 

• 575 K 

• 350 K 


“Test cell is cylindrical, plexiglass; dimensions given are Height x Diameter. 


D. Two-Phase, Drop Tower Tests 


Test Cell a 

Heater Size 

Heater Temp. 

Liquid 

Liquid Fill Levef 

Gas/ Vapor Tracers 

9 10 cm x 10 cm 

9 1 cm 

• 2.5 cm 

• 5 cm 

• 800 K 

• 575 K 

Silicon Oil 

9 0.5 

• 0.25 

• 0.125 

• None 

• Tobacco smoke 

• 10 cm x 10 cm 

• 1 cm 

• 2.5 cm 

• 5 cm 

• 800 K 

• 575 K 

R113 

• 0.5 

• 0.25 

• 0.125 

• None 

• Tobacco smoke 

• R113 vapor 

• 5 cm x 5 cm 

• 1 cm 

• 2.5 cm 

• 800 K 

• 575 K 

Silicon Oil 

• 0.5 

• 0.125 

• None 

• Tobacco smoke 

• 5 cm x 5 cm 

• 1 cm 

• 2.5 cm 


R113 

• 0.5 

• 0.125 

• None 

• Tobacco smoke 

• 1 cm x 1 cm 

• 5 mm 


Silicon Oil 

• 0.5 

• None 

• Tobacco smoke 

• 1 cm x 1 cm 

• 5 mm 

• 575 K 

R1 13 

• 0.5 

• None 

• Tobacco smoke 


“Test cell is cylindrical, plexiglass; dimensions given are Height x Diameter. 
b Liquid fill level is expressed as fraction of container height. 


E. Two-Phase, Space-Based Tests 


Test Celt 3, 

Heater Size 

Heater Temp. 

Liquid 

Liquid Fill Levef* 

Gas/ Vapor Tracers 

9 TBD 

• 2.5 cm 

• 800 K 

R113 

0.125 


(10 cm x 10 cm) 

• 5 cm 

• 350 K 


0.5 

(R113 vapor) 


“Test cel] is cylindrical; dimensions given are Height x Diameter. 
b Liquid fill level is expressed as fraction of container height. 
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As indicated in Table 1A, many parametric studies need to be conducted in 
the 1-g experiments. The order in which these tests should be performed is given 
in Table IB. This table specifies a large number of 1-g tests since there are several 
f ossible combinations for each test cell in the table. For example, nine tests are 
specified for the single-phase, 1-g experiments with the 10 cm x 10 cm test cell 
because three heater temperatures are specified for each of the three heater sizes. 
Even more combinations are specified for the two-phase, 1-g tests. However, some of 
these tests may be omitted from the matrix if experiments indicate that they may be 
superfluous for completing a parametric study. The following three guidelines should 
be adhered to lor each test conducted: 

1. Each test should be repeated at least twice to assure reproducibility. Success is 
measured by achieving approximately equal results (differing by no more than 
10%) two consecutive times. 

2. Some te^ts need to be done where the laser path is modified (i.e., reduce the 
laser power or bring the laser in from the bottom). 

3. Scale measurements are required prior to the two-phase tests in order to correct 
for parallax error. 

The heater temperatures 350 K and 800 K were selected to provide an order of magni- 
tude range of 0gLAT. The 575 K heater temperature was selected as an intermediate 
setting. Two heater sizes or two different heater emissivities should also be used with 
the test cells to obtain a variation of the radiation heat transfer effect. The interaction 
of droplet/particle tracers and thermocouples with the fluid flow will be evaluated in 
Ihe 1-g tests. Single-phase experiments will also be conducted in the 1-g tests to help 
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verify the computational code. The same test cells can be used for these experiments 
as will be used for the two-phase tests. 

Only three variables will be changed in the space-based experiments: Heater 
temperature, heater size, and liquid fill level. Of these, only the former is deemed 
a minimum requirement (in order to study the required range in 0gLAT). Varying 
the other two parameters is desirable. Two options for each of these three variables 
implies eight different tests planned for the Shuttle experiment. 


3.2 Measurement Requirements 


Measurement requirements vary between the three experiment settings and are 
summarized in Table 2. While it is highly desirable to obtain detailed temperature 
and velocity field measurements throughout the flow field, it is also very difficult and 
expensive to do so. Furthermore, such detail is not necessary to achieve the goal of the 
experimentation: to verify experimentally the results of the computational research. 
For this reason both “minimum” and “preferred” requirements are specified in this 
section. Most of the detailed measurements will be taken in the 1-g experiments. The 
reduced gravity tests therefore will require a minimum of flow visualization and tem- 
perature measurements to validate the simulated results. The measurement require- 
ments below were developed for the minimum and maximum heater temperatures 
which will be used in the experiments (namely 350 K and 800 K). Requirements for 
intermediate heater temperatures (such as 575 K) may be linearly interpolated. 
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Table 2: Science Requirements 


Design of experimental apparatus : 


Hem Minimum Requirement Preferred Requirement 


Range of Grashof 
number covered 

10 3 -10 6 

10 2 or 10 3 - 10 7 

Range of PgL&T covered 

3 orders of magnitude 

5 orders of magnitude 

Recommendation: Vary 

(1) G-level by different experiment settings 
to obtain 3 order of magnitude range. 

(2) Heater temperature. 

(3) Size of test cell. 

Temperature variation 
of heater 

± 1%/second at 1 
heater temperature 

Same for 2-3 heater temperatures 

Liquid fill level 

1 setting at 50% of 
cell height 

2-3 settings to study effect of 
cr in shallow versus deep pools 

Size and/or emissivity 
of heater 

1 setting 

2-3 settings to vary radiation 
heat transfer effect 

Measurement requirements: 



Item 

Minimum Requirement 

Preferred Requirement 

Velocity 

Qualitative full field 
flow visualization. 

Film synchronized to 
0.5 second with 
the run time of 
the experiment. 

Quantified flow velocities; 

• With 800 K heater a : 

1-10 cm/sec ±1-2 mm/sec 

• With 350 K heater 0 : 

0. 5-5.0 mm/sec ±0. 1-0.2 mm/sec 
See Section 3.2, requirement #1 
for spatial resolution 


requirements. 


Film synchronized to 
0.005 second with 
the measured data. 

“Note: Anticipated velocities in liquid phase more than 2 mm below surface are on order of 0-0.1 mm/sec 
(may be too difficult to resolve quantitatively) for experiment times of less than 10 minutes. 
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Table 2: Science Requirements (continued) 


Measurement requirements: (continued) 


Item 

Minimum Requirement 

Preferred Requirement 

Temperature 

Heater, ambient, and wall 
temperatures; also fluid 
temperatures for 1-g tests. 

• With 800 K heater: 
280-810 K ±0.5°C 

• With 350 K heater: 
280-355 K dtO.rC 

Data rate: 1 Hz 

In addition, fluid temperatures 
especially along liquid surface. 
See Section 3.2, requirement #2 
for spatial resolution 
requirements. 

Pressure 

Ambient and gas phase 
in test cell 11.7-17.7 
psia ±0.2 psia 
Data rate: 1 Hz 

Same 

Gravitational 

acceleration 

In space-based 
experiment for normal 
gravity direction 
10 -3 -10 -4 ± 5 x 10- 5 g 
Data rate: 5 Hz 

All 3 principal directions 
for space-based experiment 
plus normal gravity direction 
for 1 drop tower run 
(latter 10* 5 ± 10" 6 g) 

Data rate: 10-100 Hz 

Interfacial position 

Measure minimum and 
maximum height of 
liquid surface to ±4% 
of radius (half-length) 
of test cell. 

Measure position and 
shape to ±4% of radius 
(half-length) of test cell 
every 1 sec until change 
is less than 1 mm/sec. 

Experimental conditions requirements: 


Item 

Minimum Requirement 

Preferred Requirement 

Ambient temperature 

Shall be between 240 K 
and 320 K 

300 K ±5 K 

Ambient pressure 

14.7 ±3 psia 

Same 

G-jitter 

Shall not exceed 10“ 4 g 

Shall not exceed 5% of the 
normal gravitational 
acceleration. 

Interfacial position 

Local angle of liquid 
surface shall not 
exceed 0.04 radian. 

Local angle of liquid 
surface shall not 
exceed 0.02 radian. 
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1. Velocity Measurements 


Minimum Requirement : Flow visualization shall be used to measure qual- 
itatively the flow pattern with a full field of view in each of the three experiment 
settings. Flow visualization is required in one (1) vertical plane of the test cell 
parallel to the expected two-dimensional or axisymmetric motion. The flow 
visualization record should be synchronized to approximately 0.5 second with 
the run time of the experiment. 

Preferred Requirement : In addition to the above, it is desirable to quantify 
the velocity measurements. Experiments using the 800 K heater should measure 
velocities on the order of 1-10 cm/sec within an accuracy of ±1-2 mm/sec. 
Experiments using the 350 K heater should measure velocities on the order of 
0. 5-5.0 mm/sec within an accuracy of ±0. 1-0.2 mm/sec. Anticipated velocities 
more than 2 mm below the liquid surface are on the order of 0-0.1 mm/sec. 
Therefore it may be too difficult to resolve quantitatively the velocities in this 
regime. Film from the flow visualization should by synchronized to 0.005 second 
with the measured data. The fluid velocity measurements should be spatially 
resolved with the following accuracy: 

Near inside of wall (gas phase): 

±1 mm within 5 mm of wall for 1-g experiments. 

±1 mm within 10 mm of wall for reduced gravity experiments. 

Near liquid-vapor interface: 

±2 mm within 5 mm of interface (in gas phase) for 1-g experiments. 

±2 mm within 10 mm of interface (in gas phase) for reduced gravity experi- 
ments. 

±2 mm within 2 mm of interface (in liquid phase) for all experiments. 
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Elsewhere: 

±4 mm (gas phase). 

±10 mm (liquid phase). 

Quantitative flow measurements should be taken at a frequency of 1 Hz through- 
out each experiment. Data should be obtained for the full field of flow if possible 
(measurements along the liquid surface are the most desirable). At some point 
during each test it is required to check for axisymmetry. 

Rationale : Qualitative flow visualization will reveal the number and shape 
of the recirculation cells. Such information is necessary to compare with the 
predictions of the computational code as well as with previous test results from 
related experiments. Quantifying the flow field velocities is highly desirable 
in order to further validate the model. Without quantified flow velocities, a 
detailed comparison of the dominating mechanisms cannot be made. The ex- 
pected orders of magnitude of the flow velocity for the different tests are based 
on computer simulation results. The desired accuracy of the measurements 
is based on a deemed allowable ±10-20% error in the measurements. This is 
somewhat arbitrary. 

Spatial resolution requirements were derived from anticipated velocity gra- 
dients based on simulated results. Where the gradients are large, a higher spatial 
resolution is required. The accuracies cited in the preferred requirement are on 
the order of a grid size used for the respective regions in simulations run on the 
CRAY computer. If the true velocity gradient is known, then the error in a ve- 
locity measurement due to an error in spatial resolution (say in the y-direction) 


This equation shows that for an allowable error in velocity measurement due 
only to an error in spatial resolution, the required accuracy on the location of the 
measurement is inversely proportional to the velocity gradient. The quantified 
requirements given above yield allowable errors in the velocity measurements 
due to uncertain spatial resolutions in all regions except those very close to 
a wall (where the no slip condition results in very large velocity gradients). 
Large velocity measurement errors will exist for these regions as more stringent 
requirements than ±1 mm do not seem practical. The preferred film synchro- 
nization requirement reflects a ±0.5 mm contribution to the error in spatial 
resolution measurements for 10 cm/sec velocities. 

2. Temperature Measurements 

Minimum Requirement : In all experiments, temperature measurements 
of the heater surface and the test cell walls are required. The temperature of 
the ambient environment shall also be measured in the 1-g and space-based 
experiments. The locations of the heater surface and wall temperature mea- 
surements are shown on Figures 3.1 and 3.2. In these figures, the locations of 
the thermocouples are specified more as a guideline than as a precise require- 
ment. The most important temperature measurements are the ones adjacent 
to the liquid surface and also those on and adjacent to the heater. The other 
thermocouples which are shown on the figure are either equally spaced along 
the walls of the gas phase in one-half of one plane, or are situated in order to 
check for axisymmetrv of the temperature field. Although not shown in the 
figures, future tests may also incorporate a rake with thermocouples mounted 
along the centerline immediately below the liquid surface (similar to the design 
used in Reference [13]). The accuracy of the temperature measurements shall 
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|*-(R-RH)/3 » | ^ RH »| 



Notes : (t ) Thermocouples are on the inner wall at r or y locations specified in above 
figure except where noted below. 

(2) The two thermocouples at r « RH/3 and r = 2RH/3 are on top of the heater. 

(3) The two thermocouples marked with the symbol "A" are on the outside 
of the test cell wall. These thermocouples are not required for drop 
tower tests. 


(4) HG = height of gas phase, HL = height of liquid phase, R = inner radius 
of test cell, RH = radius of heater assembly. 


Figure 3.1: Thermocouple locations for heater surface and wall temperature measure- 
ments (cylindrical geometry). 
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Notes: (1) Thermocouples are on the inner wall at x or y locations specified in above 
figure except where noted below. 


(2) The two thermocouples at x ■ LH/3 and x - 2LH/3 are on top of the heater. 


(3) The two thermocouples marked with the symbol "A" are on the outside 
of the test cell wall. These thermocouples are not required for drop 
tower tests. 

(4) HG = height of gas phase, HL - height of liquid phase, 2L = inner length 
of test cell, 2LH = length of heater assembly, W - width of test cell. 


Figure 3.2: Thermocouple locations for heater surface and wall temperature measure- 
ments (planar geometry). 
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be at least ±0.5°C for tests with the 800 K heater, at least ±0.3°C for tests 
with the 575 K heater, and at least ±0.1 °C for tests with the 350 K heater. 
The wall temperature measurements shall be resolved spatially to an accuracy 
of ±1 mm. The frequency of the measurements shall be at least 1 Hz. 

Preferred Requirement : In addition to the above, it would be desirable 
to obtain fluid temperature measurements, particularly near the liquid surface. 
The fluid temperature measurements should be resolved spatially to the follow- 
ing accuracy: 

Below top of cell: 

±1 mm within 5 mm of wall for 1-g experiments. 

±1 mm within 10 mm of wall for reduced gravity experiments. 

Near side walls (gas phase), except where specified above: 

±2 mm within 5 mm of wall for 1-g experiments. 

±2 mm within 2 mm of wall for reduced gravity experiments. 

Near liquid- vapor interface: 

±2 mm within 5 mm of interface (in gas phase) for 1-g experiments. 

±2 mm within 10 mm of interface (in gas phase) for reduced gravity experi- 
ments. 

±2 mm within 2 mm of interface (in liquid phase) for all experiments. 

Elsewhere: 

±4 mm (gas phase). 

±10 mm (liquid phase). 

Rationale : Again, quantitative measurements are required to validate the 
computational model. Such information would supplement the quantified flow 
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velocity measurements. The temperature of the heater surface must be mon- 
itored and recorded to provide input to the computer code. The required ac- 
curacy of the temperature measurements is based on the expected temperature 
rise along the inner wall of the test cell in the gas phase and along the liquid 
surface. Simulations indicate that these temperatures may rise only about 1°C 
over the course of one hour with the 350 K heater. Hence for this case a mini- 
mum allowable precision is ±0.1 °C. With the 800 K heater the inner wall may 
increase a few degrees before a quasi-steady state condition is reached. Hence 
±0.5°C would be approximately ±20% of this temperature rise. The required 
accuracy for temperature measurements with the 575 K heater was interpolated 
from the above requirements for the 350 K and 800 K heaters. The minimum 
time between temperature measurements (1 s) is commensurate with the time 
required for the inner surface temperature of a plexiglass wall to increase PC 
in response to an 800 K 2.5 cm diameter heater. This is also of the same order 
as the time required for thermal diffusion between two node points in the con- 
tainer wall in the computational model (based on 10 node points in a 1/4” thick 
glass wall). Similar to the velocity requirements, spatial resolution requirements 
were derived from anticipated temperature gradients based on simulated results. 
Very close to the heater large temperature measurement errors will exist due 
to the large temperature gradients. For example, simulations indicate that gra- 
dients on the order of 100-200°C/mm will exist 1-2 mm below the heater. In 
this region an error of ±1 mm in spatial resolution whould therefore result in 
temperature measurement uncertainties of 100-200°C. Such a large uncertainty 
in this region may have to be accepted, as more stringent requirements than 
±1 mm do not seem practical. 
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3. Pressure Measurements 


Minimum Requirement : It is required to measure the ambient pressure 
and the average pressure of the gas phase in the test cell. The accuracy of the 
measurements shall be ±0.2 psia. The frequency of the measurements shall be 
at least 1 Hz. 

Rationale : Ambient pressure measurements are useful in assessing the 
possibility that gas leakage occurred during an experiment due to a pressure 
difference across a small space between the heater and the container or between 
any other small crevice of the test cell. The average pressure inside the test 
cell and the bulk temperature of the gas is required to determine the total 
vaporization rate of the liquid as a function of time. The required accuracy of 
±0.2 psia represents approximately a 10-15% change in temperature of the gas 
inside the container or a 10-15% change in ambient pressure. 

4. Gravitational Acceleration Measurements 

Minimum Requirement : Gravitational acceleration levels shall be mea- 
sured in the space-based experiments. As a minimum, the acceleration in the 
plane of the normal gravitational force is required. The accuracy of the mea- 
surements shall be ±5 x 10 -5 g (approximately ±5 x 10 -4 m/s 2 ). Data should 
be collected as frequently as practical (data rate shall be at least 5 Hz). 

Preferred Requirement : Gravitational acceleration levels shall be mea- 
sured in the three principal directions of the space-based experiments. The accu- 
racy of the measurements shall be ±5 x 10 -5 g (approximately ±5 x 10~ 4 m/s 2 ). 
In addition, the acceleration should be measured in the plane of the normal 
gravitational force during one drop test run. The accuracy of this measurement 
shall be ±10~ 6 g (approximately ±10 -5 m/s 2 ). The data collection rate for the 
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gravitational acceleration shall be at least 5 Hz for the space- based experiments, 
and at least 10 Hz for the drop tower test. 

Rationale : The computational code can incorporate an oscillating com- 
ponent of the normal gravitational force in order to simulate g-jitter which is 
encountered on manned space vehicles. These random fluctuations occur at 
high frequency, so the data collection rate must be correspondingly high. A 
minimum rate of 5 Hz (data every 0.2 second) was selected because this is 
comparable to the time step used in the simulation for a low-gravity run after 
achieving a quasi-steady state. A greater frequency of 10-100 Hz would be 
more desirable to represent more accurately the fluctuating g-jitter and to cor- 
respond with the smaller time steps required in the simulation prior to achieving 
a quasi -steady state. Gravitational acceleration levels orthogonal to the normal 
gravitational force would be desirable to compare with the magnitude of the 
normal gravitational acceleration. Such measurements are necessary to assess 
the effect of these acceleration components on the flow. The data collection 
rate for these orthogonal gravitational acceleration components need not be as 
high (about 1 Hz) as the normal gravitational acceleration rate. The accuracy 
of the acceleration measurements represent ±10% of the expected mean normal 
gravitational acceleration levels. 

5. Measurement of Interfacial Position 

Minimum Requirement : A qualitative picture shall be taken of the portion 
of the liquid surface which meets the front side wall (facing the camera). The 
minimum and maximum height of the liquid surface must be measured with an 
accuracy of ±4% of the radius (half-length) of the test cell. 
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Preferred Requirement : The position and shape of the liquid surface shall 
be measured with an accuracy of ±4% of the radius (half-length) of the test 
cell. If an initial change in the position of the liquid surface occurs (as expected 
in drop tower tests), then the position and shape of the liquid surface shall be 
measured every 1 second until the position of the interface changes by less than 
1 mm over a period of 1 second (or until the drop is complete for the case of 
the drop tower tests). 

Rationale : The height of the liquid pool must be determined. The required 
accuracy of ±4% of the radius (half-length) of the test cell is comparable to the 
related requirement of allowable liquid surface deformation. It is also desirable 
to obtain the quantitative shape of the liquid surface so that thermocouples may 
be placed accurately right above the liquid surface to obtain important temper- 
ature measurements along the interface. If necessary, quantitative information 
on the shape of the liquid surface may also be used in a future modification of 
the computational model which allows for a curved gas/liquid interface. Since 
obtaining quantitative information on the shape of the liquid surface may be 
exceptionally difficult to accomplish, this has been cited as a “preferred” rather 
than a “minimum” requirement. 


3.3 Experimental Conditions 


Some external conditions must be controlled to minimize their effects on the 
experiment. The requirements for the external conditions are given below. In general 
these requirements are not very restrictive because the computational model incorpo- 
rates environmental conditions such as ambient temperature and pressure and normal 
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(vertical) gravitational force. Hence measurements of these parameters is the only 
strict requirement. However, as noted below, operating the experiments under “nor- 
mal” external conditions is the preferred requirement. 

1. Ambient Temperature 

Minimum Requirement : The ambient temperature shall not exceed 320 K 
nor be less than 240 K during the experiment. 

Preferred Requirement : The ambient temperature should be maintained 
at 300 ±5 K throughout the experiment. 

Rationale : Since liquid R 1 1 3 vaporizes at 320.7 K, a maximum ambient 
temperature of 320 K was selected. The lower bound of 240 K was selected 
because little thermophysical property data exists for liquid R.113 below this 
temperature. No restrictions on the shape of the temperature profile were 
made since ambient temperature measurements can be incorporated into the 
computational model. A preferred requirement of maintaining the ambient 
temperature at 300 ±5 K was selected because simulations show that the liquid 
surface and wall temperatures of the experiment change by not more than ±1% 
over 1 hour in the case of a 305 K ambient temperature versus a 300 K ambient 
temperature. 

2. Ambient Pressure 

Minimum Requirement : The ambient pressure shall be maintained at 
14.7 ±3 psia throughout the experiment. 

Rationale : To minimize leakage of R113 vapor through imperfect seals in 
the container, the ambient pressure should be maintained within about 20% 
of the initial gas pressure in the test cell. Pressure fluctuations should not 
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occur over a period of time less than that required for the gas pressure to 
equilibrate. This minimum fluctuation time is TBD. Another criterion for the 
ambient pressure requirement is based on keeping the flow velocity through a 
crack (orifice due to imperfect seal) in the test cell much less than the flow 
velocity in the test chamber. This analysis is TBD. 

3. G-jitter 

Minimum Requirement : The g-jitter level shall be kept at 10 -4 g or less 
during the test in order to keep the amount of free surface deformation below 
1.5 mm. 

Preferred Requirement : The g-jitter level should be kept below 5% of the 
normal gravitational acceleration during the experiment. 

Rationale : The assumed mean gravitational level of 5 x 10~ 4 g on the 
Shuttle is acceptable for the space-based experiments in order to achieve the 
required range of (3gLA.T and the Grashof number. The effect of g-jitter during 
the space-based experiments not only changes the buoyancy force in the exper- 
iment, but also affects the shape of the liquid free surface. By measuring the 
gravitational force during the experiment, the input of gravitational force in the 
computer simulation can be changed accordingly. However, the computational 
code still assumes a flat liquid-gas interface. Deformation of the liquid free sur- 
face will be minimized by the liquid-gas interface control system (the internal 
flanges — see Appendix A). However, the g-jitter level should be kept at 10 -4 g 
or less at any time during the test to keep the local angle of the liquid surface 
below 0.04 radian. The preferred requirement of a 5% g-jitter level is desired in 
order to minimize the effect of g-jitter and thus simplify the analysis. 
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4. Interfacial Position 


Minimum Requirement : The local angle of the liquid surface with respect 
to the plane that is perpendicular to the side walls shall be no greater than 
0.04 radian. 

Preferred Requirement : The above local angle of the liquid surface shall 
be no greater than 0.02 radian. 

Rationale : Based on static considerations, the liquid surface deformation 
must be minimized so that the vertical component of liquid surface velocity is 
negligible compared to the horizontal component. The local angle of the liquid 
surface is 

e “ Une= v: = Wr-6 

where L is the radius (half-length) of the test cell, R is the radius of curvature 
of the liquid surface, and 8 is the liquid surface deformation with respect to 
the height of the liquid in contact with the side walls. It can be shown that 

if R » L, then 

8 

0 < 1 =* - < 1 

Therefore the requirements above on the maximum local angle of the liquid 
surface represent approximately a 4% and 2% maximum allowable liquid surface 
deformation with respect to the radius (or half-length) of the test cell. 
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3.4 Preliminary Hardware Concepts 


The ground-based and space-based experiments consist of the following hard- 
ware: 

1. Test cell: Plexiglass. 

2. Fluids (liquid and gas): Silicon oil, R113, or water (liquid); air (gas). 

3. Heater: Nichrome wire in ceramic cement. 

4. Flow visualization system: Laser light sheet with tracer particles. 

5. Temperature measurement system: Thermocouples (microgravity and 1-g), in- 
terferometry (1-g only). 

6. Liquid-gas interface control system: Internal flanges. 

7. Data acquisition and control system. 

Hardware concepts for the ground-based and space-based experiments are given in 
Appendix A. 
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Chapter 4 

Experimental Timeline 


4.1 1-g Tests 

No automation is required in the 1-g tests. A typical sequence of events is as 
follows: 

1. Fill test cell to desired level with liquid and tracer particles. 

2. Preheat heater until it reaches a steady temperature. This may be done either 
away from the test cell or in place. 

3. Install heater and entire top of test cell. 

4. Turn on laser for flow visualization system. 

5. Turn on camera(s) for recording system. 

6. Begin data collection. 

7. Use second light sheet halfway through experiment for 2 minutes to check for 
axisymmetry (see Section 3.2). Replace light sheet into original position after 
this time. 
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8. End test run when surface temperature has readied a quasi-steady state (see 
note below) or after 1 hour. 

As discussed in [1], a quasi-steady state condition is characterized by stabilized 
thermal and dynamic boundary layers. This occurs when the temperature and veloc- 
ity on the liquid surface vary very slowly. Typically this takes between 1-10 minutes. 
A truly steady state condition would require several hours to reach and is not nec- 
essary for the purpose of the experiment. Therefore a run time of 1 hour has been 
selected. It should be noted that smoke tracers or tracer droplets remain in suspension 
for only approximately 10 minutes. Since injecting more vapor tracers will interfere 
with the flow, the 1 hour runs will only have vapor tracers for the first 10 minutes. 


4.2 Drop Tower Tests 

The procedure for the drop tower tests is slightly different than for the 1-g tests 
because the experimental package hangs above the drop tower for at least 3 minutes 
prior to the drop. During this hang-time the experiment is enclosed in a capsule and 
cannot be touched. Therefore an automated system must be used to inject tracer 
droplets into the gas phase immediately before the drop. A remote means is also 
required for turning on the heater while it is hanging in the drop capsule. The 
sequence of events is as follows: 

1. Install heater and entire top of test cell. The heater should be turned off and 
at ambient temperature at this point. 

2. Fill test cell to desired level with liquid and tracer particles. 
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3. Install experimental package into drop capsule. Raise capsule to its hanging 
position. 

4. Turn on camera and monitor liquid in test cell. 

5. After liquid becomes quiescent, turn on heater (remotely) and begin data col- 
lection. 

6. Monitor heater temperature. Drop the experiment when the heater reaches the 
desired temperature. 

7. Data collection ends after approximately 2-3 or 5-6 seconds (depending on 
which drop tower facility is used) when drop is complete. 

4.3 Space-Based Tests 

The procedure for these tests is the same as for the 1-g tests, except that the 
liquid fill and tracer droplet/particle injection and stirring are done by automated 
systems. A crewmember may be required only to turn on the heater, monitor its 
temperature until it reaches a quasi-steady state, turn on and monitor the liquid 
fill system until the desired liquid fill level is achieved, and then press one button 
to simultaneously: inject and stir tracers, activate light sheet, and turn on data 
acquisition system. Each experiment will take 1 hour after the initial 20 minute 
set-up time. 
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Chapter 5 

Planned Data Analysis 


The following data collected in the experiments will be used for analysis: 

1. Video or movie films as the flow observation record. 

2. A tape (memory size TBD for the three experiment settings) with the recording 
of the thermocouple outputs. 

3. G-jitter and acceleration data during the tests (space-based only). 

4. Comments of the test operators. 

The films will be reviewed to observe general flow patterns and free surface behaviors. 
The particle motions will be studied by a motion analyzer to quantify the flow visual- 
ization at various times. The thermocouple outputs will be converted to temperatures. 
The temperatures of the heater, the ambient air, and the fluid at various locations 
(ground-based tests only) will be plotted as a function of time. Wall temperatures 
at various locations will also be plotted against time. For the 1-g tests, the surface 
temperature distribution will be plotted for various times. Results of the data analy- 
sis will be compared with the predictions of the computational model. Ground-based 
test results will be used to aid in the design of the space-based experiments. 
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Chapter 6 

The Numerical Model 


6.1 Description of Numerical Model 

The finite-difference procedure applied to the governing equations in the gas 
and liquid phases utilizes the SIMPLE algorithm [3] with the SIMPLEC modification 
[14]. The important operations of this algorithm, in the order of their execution, are 
(taken from [3]): 

1. Guess the pressure field p*. 

2. Solve the momentum equations to obtain u* and v‘ . 

3. Solve the p' equation. 

4. Calculate p by adding p' to p*. 

5. Calculate u and v from their starred values using the velocity-correction formu- 
las. 

6. Solve the discretization equation for other <j >' s (such as enthalpy and mass frac- 
tion of species) if they influence the flow field through fluid properties, source 
terms, etc. 
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7. Treat the corrected pressure p as a new guessed pressure p', return to step 2, 
and repeat the whole procedure until a converged solution is obtained. 

The numerical model uses primitive variables (u,t>,p, h), a staggered mesh [3], and 
variable thermophysical properties (density, viscosity, etc.) in both phases. The 
effects of surface tension, evaporation, radiation, and wall thermal inertia are incor- 
porated into the computational code. The user has the option of using an explicit 
one-dimensional or two-dimensional wall conductivity calculation (the latter takes 
into account heat transfer between the outer walls and the environment via Newton 
Cooling). All surfaces are assumed to have an emissivity of 1.0. The calculation of the 
radiative heat transfer between elementary rings of the enclosure and liquid surface 
has been described previously by Abramzon, Edwards, and Sirignano [1], The pro- 
gram also accounts for the total pressure variations in the gas phase. This option is 
very important for the enclosed container simulations because the total gas pressure 
is continuously increasing during the heating and vaporization processes. The code 
also allows for local refining of the finite-difference grid at regions where the field 
gradients are large and better resolution of the problem is needed. 

The following variables are allowed to be changed parametrically: 

1. Liquid fill level. 

2. Container dimensions, including wall thickness. 

3. Thermal driving potential (heater temperature, position, and size). 

4. Fluid properties. 

5. Gravity level. 
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The gas and liquid phase solutions are time-split in that the solution at time t rep- 
resents the gas phase at time t and the liquid phase at time t -f |a t (where A< is 
the time step). At the liquid- vapor interface, values of liquid temperature and veloc- 
ity are used for the gas phase solution while values of gradients of gas temperatures 
and velocities are used for the liquid phase solution. Because the prediction of the 
liquid-vapor interfacial shape is complex, the model presently considers only a flat 
liquid-vapor interface. The sequential calculation procedure for each time step is as 
follows: 

• Determine wall boundary temperatures using explicit scheme: 

1. Calculate total heat fluxes (radiative + conductive) to the walls in the gas 
phase based on temperatures at time t. 

2. Calculate conductive heat fluxes to the walls in the liquid phase based on 
temperatures at time t. 

3. Use explicit scheme to calculate wall temperatures based on above heat 
flux data. Newton’s law of cooling is used as a boundary condition for 
the outside of the walls to take into account natural convection to the 
environment. 

• Determine gas phase solution at time t + At: 

1. Calculate gas phase properties at beginning of time step. Only values of 
c p are updated during the internal iteration cycle. This helps stabilize the 
solution. Values of c p are updated in order to update values of enthalpy. 
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2. Perform internal iteration loop for gas phase (using SIMPLEC algorithm). 
Boundary values of wall temperatures (from explicit solution) and liquid 
temperatures and velocities (from time t 4- ^A t) are used. 

• Determine liquid phase solution at time t + -At: 

1 . Calculate liquid phase properties at beginning of time step. Only values 
of c p are updated during the internal iteration cycle. 

2. Perform internal iteration loop for liquid phase (SIMPLEC algorithm). 
Boundary values of wall temperatures (from explicit solution) and gradi- 
ents of gas temperatures and velocities (from time t + At) are used (e.g., 
heat and stress balances at interface). The bottom of the liquid phase is 
assumed to be isothermal at the reference temperature. 

As previously mentioned, the resulting solution for time t + At actually represents the 
gas phase at time t+At and the liquid phase at time t+f At. The calculation procedure 
for the single-phase simulation is similar to the above except that a relatively thick 
wall is used at the bottom of the gas phase with the isothermal boundary condition 
being applied at the bottom of this wall. 

Calculations are typically performed with a 22 x 22 mesh in each phase for as- 
pect ratios equal to 1.0. Since two boundary nodes are used in each direction, there 
are 20 x 20 internal cells in this mesh. Denser meshes were employed to test grid sen- 
sitivity. For one simulation, three runs were made with mesh sizes of 22 x 22, 32 x 32, 
and 42 x 42. The resulting stream function contour plots after t = 100 seconds were 
qualitatively identical. Quantitatively, the local maxima of the stream function dif- 
fered by only 1% between the 32 x 32 and 42 x 42 meshes, while these values differed 
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by 8-10% between the 22 x 22 and 32 x 32 meshes. The time step was typically about 
0.05 second and required about 8 seconds of CPU time on a VAX- 11/780 computer 
for the two-phase calculations, and about 4 seconds of CPU time for the single-phase 
calculations. For the single-phase simulations presented in Chapter 7, a 22 x 42 grid 
was used with a 0.01 second time step. 


6.2 Governing Equations 

Aside from the continuity equation, which we solve for by iterating with the 
pressure correction equation, the following conservation equations are solved in nondi- 
mensional form: momentum equations (u- and v-), the energy equation (for excess 
enthalpy), and species equations (for mass fractions). Each equation is written in a 
general form. 


6.2.1 The General Conservation Equation 


The general form of the dimensional conservation equation is as follows: 


( 6 . 1 ) 


where 


<j>' = the quantity of interest (e.g., u , Y- , h\) 

r; = dimensional diffusion coefficient for <f> 
(e.g., //, k'f%) 
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dimensional source term for <f> 

r 

x' for cylindrical coordinates 

< 

1 for rectangular coordinates 


To make equation (6.1) nondimensional, use: 


t = t'/t„ <f> = x = x'/L., y = yV L., 


14 = WJU* an d Vy = VJ/t/. where U. = 

17 = r;/r. (where T. = for all <£), g = g't./U *, 

S, = /> s P'/(AP)„ t = T = TIT., and 

p.4>. p.U. 

T, = T'JT. where T[ = (A T)' = T'-T . and (AT), = T. 
By multiplying equation (6.1) by 

t. _ L. 
p.4>. p.U.<t>. 


and noting that 


we get the following 


r. 


t dll = MR* 
U.L. ~ 1/R 


P.U.L. 

nondimensional general conservation equation: 


W) 18. 

dt r dx 


pVA - -E^r* 


_i_ 

Re. 


d£ 

dx 


+ a?r^‘ 


df 


Re. 


r *tr s * 


6.2.2 Continuity Equation 


The continuity equation in nondimensional form is 



In cylindrical coordinates with d/dd = 0, this becomes 

dp 13 d . . 

Si + -rTr {prV ' ] + Ty W = ° 

The rectangular coordinate equation and the cylindrical coordinate equation can be 


combined by 

dp 1 d d . 

Si + vTx (rpu) + S^ {fm) = 0 (64) 

where r = x for the cylindrical geometry, r = 1 for the planar geometry, u = V r = V x , 
and v = V y . This designation will be assumed for the remainder of this document. 
Note also that equation (6.4) is in the general form of equation (6.3) with 


<j> = 1, = 0, and S 4 = 0 


6.2.3 Momentum Equations 


Omitting the single quote (') superscripts, the dimensional momentum equation 


is 


DV dV - -. 

P ~Dt = P ~di + = - VP + t V ^ + P9 


(6.5) 


where [V*f] is not a simple divergence because of the tensorial nature of f; t is the 
stress tensor defined for the cylindrical geometry as (convention opposite to [15]): 

( 6 . 6 ) 


2*1 - li}-9 

dx 3 


T rr — T xx — p 

Tee = p 
(note Tge = 0 for planar geometry) 


u 2 - - 
2- _ f V-V 
L r 3 


Tyy — P 


dv 2 ~ .-i 

i% ~ i 


TtO T$ r T$y — Ty0 — 0 

T yx ~ Try ~ ft 


dv du 
dx dy 


(6.7) 


(6.8) 

(6.9) 

( 6 . 10 ) 
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The nondimensional form of the momentum equations will now be treated separately 
for the r (or x) and y directions. 


u-momentum equation 


Equations (6.5) and (6.6)-(6.10) for the r (or x) direction give 


du du 


du 


dt +U fc +V d^)- 


dP 


dx 


+ 


1 ^ < \ T 68 
;^ (rr “ ) - t + sf 


=0 

+ P 'sT 


( 6 , 11 ) 


The left hand side of this equation can be re-written with the aid of the continuity 
equation as 


=o 


f du du du N 

7T7 + u ^ h V ~^T 

ut ux Oy i 


+ 


dp Id. d . 


U = 


djpu) 

at 


+ u 


1 d . , d , ' 

-rTx {rPU) + Ty (pV) 


du 


du 


+ pU d~x + pV Ty 


djpu) 

dt 


+ 


u d 
r dx 


(rpu) + pu 


du 

dx 


+ 


u—(pv) + pv 


du 

dy. 


djpu) 

dt 


m l— 

^ r dx 


( rpu • u) + 


d_ 

dy 


(pvu) 


djpfl 

dt 


, 1 d , d . 

+ ~—(rpu4) + ^ W) 


( 6 . 12 ) 


where 4> = u. 


Substituting equation (6.12) into equation (6.11) gives 


d(p<f>) Id. ^ d . dP 

_ + __ (rW) + _ ( ^ ) = __ + 


IiL( rT ) 

rdx { xx) 


To get this equation into the standard form, we subtract 


i d 


du 


du 


_ ( r P a_ I a.. ( P a.. 


r dx 


dx 


dy V 


1~88 dtyx^ 

r dy 
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from both sides of the equation to get 


d(p<t>) 1 9 

dt r dx 


r (put - ^)] + | (pvt - p^j = 


_dP_ l_d_ 
dx + r dx 


r r T 


du\l d ( du\ Tgg 

dx ) dy \ Tvx ^ dy) r 

Denoting the right hand side of this equation as S ^ and letting 


<f> = u, T = y 


the above dimensional equation can be multiplied by 

U = L. 

p.4>» p.u*<t>. 


to get the nondimensional u-momentum equation: 


dt 


1 d_ 

^ r dx 


pu<f> — 


Re , dx . 


d_ 

dy 


i^pv<i> — 


_ s'+t. 

Re, dy ) p,<j > . 


(6.13) 


Therefore the nondimensional n-momentum equation is in the general form of equa- 
tion (6.3) with 



(6.14) 


where the subscript * denotes a reference quantity and the superscript ' denotes a 


dimensional quantity. The source term 5^ in equation (6.14) will be defined in terms 
of nondimensional quantities in Section 6.3. 
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r-momentum equation 


Equations (6.5)— (6. 10) for the y-direction give 


P 


( dv dv dv\ 

\m+ u Tx + v T y ) 


dP 

dy 


+ 


\jL( rr \ i d T w' 

r&t ( iv) + dy 


+ P9y 


(6.15) 


where all quantities above are dimensional. Again, the left hand side of this equation 


can be rewritten using the continuity equation as 


Dv 


P 


Dt 


djpv) 

dt 

d(pv) 


+ 


vd dv 

-- (rpu)+l>u - 


Id. . d . 

at + - rx (rpu ■ ») + Ty (pv 

8(p*) 18, ^ 8 

ST + + 


d . dv 

v Y ipv)+,m Y. 

v) 


Hence it can easily be seen that the nondimensional u-momentum equation can be 
written in the general form of equation (6.3) with 


<f> = v, = /x, and 


= — {■ 

p.u. { 


dP' 1 d 

T + ~ 


dy' 


+ W [T 


r dx' 
.dv 


r \ T ry-P 


,dv_ 

dx 


yy 


l ‘W i+ps 


I / 

y 


(6.16) 


Note that in equation (6.16), g' y = — \g'\ for the given orientation (e.g., for earth 
gravity g' y = -9.81 m/s 2 ). 


6.2.4 Energy Equation 


The thermal energy equation is solved in terms of the excess enthalpy of the 
mixture. The derivation of this equation appears in Appendix B while the final result 
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[in the form of equation (6.3)] is given below. The excess enthalpy of the species 
(h{ .excess) an< ^ °f ^he mixture (h\) is defined by 


fc i, excess 


= c;,(r'-r') = c; t; 


and 


where 


h\ = Y,Y'A.i = cX = h'-h°' 

X 

= h i + /i, excess = "h 
/if = h t (T° ) = formation enthalpy 
T° = T. = 300 K (reference temperature) 
^ 

^ = EW 

t 

(rrw)£ c '» {T ' )dT ' 


(6.17) 


(6.18) 


C Pl 


C » = 


The dimensional energy equation written in the standard form of equation (6.3) using 
the mass-weighted average velocity components ( u and v) is 


d . 1 d 

W i/,h ' ) + 7d7 


i (, i i k' dh\ \ 

• ('“*■- if er) j 


where 


d / , , t , k' dh[\ 

+ wr h '-%w mS * 

s ‘* = w + + *'• {"'W' [bw> - ^ v?;] } 

+,*■ + & -\e£(ZY;VrV;) 

-V'- {/ + |r;K(K.V7')]| + v- 5 p'i? 


(6.19) 


(6.20) 
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and 




Le = 


V 

*} 

k'/c, 


p _ 


Qv — 


p'% 

dQ' 


dV' 


heat released by chemical reaction (6.21) 

(calculated at temperature T° ) 

binary diffusion coefficient 

Lewis number of mixture (6.22) 

viscous dissipation function (see [16], eqn A. 16) 
volumetric heat source term 

(energy production per unit volume per unit time) 


f/' Y i 

* - Y ; 


(6.23) 

v'a = E y& 

(6.24) 

v { = v' A + v; 

(6.25) 


Note that v { is the local average velocity of species i (relative to stationary coordinate 
axes), v a is the mass-weighted average velocity of the fluid mixture, and Vi is the 
diffusion velocity of species i. The nondimensional form of equation (6.19) is 


djpf} ld_ 

dt r dx 


JL( A r * Q£\ tl a' 

+ dy\ pv4> Re.dy) p.c Pt T m S * 


with 


<f> =s h\ = 

c„,T. 


and 


r„ = 


k'/c r 


(6.26) 


t. 


The source term 

- 7T¥ S * 

in equation (6.26) [see also equation (6.20)] will be defined in terms of nondimensional 
quantities in Section 6.3. 
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6.2.5 Species Equation 


One species equation must be written for each species in the gas phase: fuel 
vapor, oxygen, and products. These equations are used in the code if vaporization 
effects are included (user option). For diffusion in a nonstationary medium, the total 
mass flux, rfi, of species i is given by (all quantities dimensional unless otherwise 
specified): 

n, = P mix Yiv’ - = pY^v' + V t ) (6.27) 

If we designate r" as the rate of production of the component (kg/(m 3 -s)], a 
two-dimensional mass balance gives 

(mass accumulated = mass produced + flux in — flux out) 

[(/>^) t+A < - W)J AxA v = 

r'” Ax Ay At + ( n x - n x+£t , x )AyAt + ( n y - n y+AV )AxAt 

Dividing by Ax Ay At and taking the limit as Ax, Ay, At — * 0 gives 

| \(pY,) = -V-n' + r" (6.28) 

Substituting equation (6.27) into (6.28) and using the general form applicable to both 
cylindrical and rectangular coordinates yields 

J \(pYi) + V'ipv'Yi - pVijVYi) = r” 
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AH quantities in equation (6.29) are dimensional. To make the equation dimensionless, 
multiply it by 


t. L. 
P.Y. ~ p.U.Y, 


to get 


d , 1 d 

Ft (pY ' ) + 7te 


r puY, - 


p. p-D’,, SY, 
p+U,L, p , dx , 


+ 


dy 


( v. ti t^ii dYi \ h. ( '"V 
^ P.U.L.' p, dy) p.Y, ( r ‘ ' 


Using the Schmidt number, Sc, = p,j p,T>„ and the nondimensionalization 


Yi = Y'/Y, (= <f> and Y, = 1) 

pV,,/Sc. = /©>. (= r*) 


we obtain the result 


d . Id 
ai W) + 'rTx 


r pu<f> — 


U_d±' 

Re, dx , 


. 9 f , T* d<f>\ __ t, / <«\' 
+ dy\ pV<t> Re,dy) ~ p.Y.V*) 

This is in the general form of equation (6.3) with 


(6.30) 


<t> = ^ 


r* = pVij/Sc, 


s* = 



(6.31) 


The computational model currently includes vaporization (and thus multicomponents 
in the gas phase) but not chemical reaction. Therefore S+ = 0 in equation (6.31). 
In the future, chemical reactions may be incorporated in the model with a simple 
one-step reaction and an Arrhenius expression for the specific rate constant. The 
source term for this one-step chemical reaction is derived in Appendix C. 
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6.3 Source Term Approximations 


The nondimensional source terms in the governing equations were defined in 
Section 6.2 in terms of dimensional variables and reference quantities. The nondi- 
mensionalization of these source terms is treated in this section. Simplified source 
terms are used in the computational model. The approximations which led to these 
simplified terms are also given in this section. 


6.3.1 Momentum Equations 


With the exception of the buoyancy term in the v - momentum equation, the 
source term of the u- and u-momentum equations are treated similarly. 


u-momentum equation 


The source term in the u-momentum equation was given by equation (6.14): 



Using the definitions of the nondimensional variables in equation (6.2), equation (6.14) 
becomes 


c _ (AP).^P , fu [l±\ r ( T _ u dn\ 

* p m Ul dx pM.L m \rdx[ \ * ^ dx ) 
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Defining 


Eu. = 


and noting that 


(An. 

pM* 

= 1/ Re. 


p.U.L. 

the nondimensional source term for the u-momentum equation becomes 


(6.32) 


(6.33) 


c _ _ r (9P\.l _l ] ± 

' * jU '{dx) + Re.\rdx 

0 ( du\ Tee' { 

f, dy) r/ 


r \t xx - ft 


du 

dx 


(6.34) 


? -momentum equation 

Before the source term in the r- momentum equation is evaluated, we must first 
examine how we define P' in the computational model (pressure splitting). The total 
pressure is divided into three parts: global pressure (P 0 ) , dynamic pressure {P[), and 
the hydrostatic pressure (p.g y y ~ —p.g y )■ Hence 

P‘ = Pl(t) + P[(r,y) - p-g'y' (6.35) 

While global pressure is used to calculate density in the gas phase, only the 
dynamic pressure appears in the momentum and pressure correction equations. In 
the u-momentum equation, the global pressure and the hydrostatic term do not ap- 
pear because only the gradient in the x-direction appears. The global pressure also 
does not appear in the r-momentum equation. However, the hydrostatic term in 
equation (6.35) contributes to the pressure gradient in the y-direction. Therefore, the 
source term in the u-momentum equation becomes [from equation (6.16)] 
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S* 


t. f 

6P [ 1 d 

t 

(r 

p.u. \ 

dy‘ + r' dx' 

T 

\ Xy * dx 'J. 



(6.36) 


In the above equation, g' y from equation (6.16) has been replaced by —g where g > 0. 
Using equations (6.2), (6.32), and (6.33), the nondimensional source term for the v- 
momentum equation becomes (note P is understood now to be the dynamic pressure) 



(6.37) 


If the user chooses to run a simulation using the Boussinesq approximation, 
then the last term in equation (6.37) is replaced by g(T — 1.0) for the gas phase 
and g/3'T.(T — 1.0) for the liquid phase (where T is evaluated at the center of the 
v- momentum control volume). 


6.3.2 Energy Equation 


The dimensional source term of the energy equation (omitting the primes on 
the right hand side of the equation) was previously given by equation (6.20) as 

s; = ^ + <?ch + V- 1 pVijT, [ifoVX) - Llvz, | 

-v- DKv;(t/-iy + jr.v'ti 7 .-' 7 .)]} + v- £ p,v, 
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This source term expression may be simplified by the following approximations: 


1. Low Mach number flow with no explosions. 

2. All species' average heat capacitances are equal and do not vary spatially for a 
given time t. 

3. Negligible viscous dissipation. 

4. No volumetric heat sources or sinks. 

5. Terms with a diffusion velocity of order two or more are negligible. 

6. The molecular weights of all species are equal. 


The first approximation above implies that 


DP 

Dt 


dP - - 
= - + (v'.V)P 


dP 

dt 


The approximation that at a given t, c Pi « constant yields 


(6.38) 


IX VK « £V(J„K) = VjfoK) = Vc p (6.39) 

t i i 

However, near the heater in the gas phase, large temperature gradients can cause 
c Pi to vary by about 5% between the adjacent cells. Thus if the approximation of 
equation (6.39) were not made, we can still introduce a different assumption by first 
writing 


XX VYi - LeVc p = 

i 

(i -Z?£(K§^vr) (6.40) 
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The derivation of equation (6.40) is as follows: 


c r=H 5 ' C P. 


Vc^B^VK+K.Vc,,.) 


and 


so that 


Vc = 

VC " dT Vi 


£c Pi vy; - LtVcp = 

i 

= Y.--,yy.-Ti n? r yy.+i:(r.^vrj 

= (i -17) Be* vy;) - (y%vr) 


From equation (6.40) it is clear that if we allow c Pl to vary for different i, we 
can still make the approximation that if Le « 1, then 


T,z r yy - £evc, = -;r (v.ffvr) 


so that 


v- {p-PijT, £(c p ,vy;) - LeV-c J J * -V- pVijTi Z ( Y % VT ) (6-41) 


The third and fourth assumptions imply that 


/x4> ss 0 


(6.42) 


and 


Qv ~ 0 


(6.43) 
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The fifth approximation that terms in which the diffusion velocity appears as 
order two or higher are neglected implies that 


1 D 
~ 2 p Dt 







(6.44 ) 


The final approximation that the molecular weights of all species are equal 
results in 

« =0 (6.46) 

U, 

since 


Using a prime (') to denote a dimensional quantity, equations (6.20), (6.38), 
and (6.41 )— (6.45) give the following nondimensional approximate source term for the 


energy equation for Le = 1: 



(6.46) 


Using 


and 


we can write 


Eu « = 


(AP). 

pM? 


Ec. — 


Ul 


1 2 


Cp.T. 


t. dP' 

P.C p , I m til 


r r dP 

Eu.tc.— 

at 


(6.47) 
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Also, since 


K = 1, P V t] = 


/ 

P v a 


and t. ~ jj- 

U * 


we can write 


±tM p ' v: ’ T ''Z { Y: ^ V r ')] = 

- 7ZJ?- l^ T ^( Y % VT ) 

pDijTi £ (i'l^vr) 


= -!-v« 

Re. 

Therefore the final nondimensional source term [equation (6.46)] may be written as 

_ _ _ dP t. , 

S* « Eu.Ec .— + z _ —Q c h 


dt p.c p .T. 


-t?\ v ^( Y % VT ) 


(6.48) 


In the current computational model, we assume that c Pi « constant, Le = 1, 
and no chemical reactions so that 

dP 

St » Eu.Ec.— (6.49) 

However, in the future when chemical reactions are included in the code, equa- 
tion (6.48) will be used as the source term of the energy equation where 

Cch^-E (KrT)' 
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Preliminary Results 


The numerical model has been tested by comparing results with a bench mark 
solution for one- phase, steady-state, buoyancy-driven flow in a square cavity [17], 
The thermal boundary conditions for this problem are differentially heated vertical 
side walls and insulated top and bottom walls. Results obtained using the SIMPLEC 
algorithm agree very well with the literature. For example, for the case of Ra == 10 6 
(Ra is the Rayleigh number). De Vahl Davis reports |— v’-’lmax = 16.772 at x = 0.154 
and y — 0.545. For this same case I find |— V’lmax = 1"*3 at x = 0.154 and y = 0.540, 
which differs from the aforementioned bench mark solution by approximately 3%. 
However, since this relatively simple test problem does not involve the effects of radia- 
tion. wall thermal inertia, or variable properties, the comparison cannot be considered 
a complete verification of the computational code. 

In order to obtain a more thorough verification of the code, it is necessary to 
compare numerical and experimental results. This involves an added complexity in 
that the computer code must accurately model the actual experiment, or vice-versa. 
If simplifications introduced in the model cannot be reproduced experimentally to 
within a prescribed accuracy, then the computational code must be modified to du- 
plicate the experimental results. As will be shown in this chapter, preliminary results 


\ 
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of a single-phase simulation agree fairly well with the corresponding experimental 
results; however some evidence exists which suggests that the current computational 
model should be modified to allow for emissivities less than 1.0. The results also 
indicate that the actual radial temperature profile of the axisymmetric heater should 
be used instead of modelling the heater as having a constant radial temperature. 

The following four single-phase, normal-gravity computer simulations have been 
performed to date for a cylindrical, plexiglass test cell (10 cm Dia. x 10 cm Height): 

1. 682 K, 1.2 cm heater with 0.4 cm thick top and side walls, and 1.2 cm thick 
bottom wall. 

2. 682 K, 1.2 cm heater with 0.4 cm thick top and side walls, and 1.2 cm thick 
bottom wall; radiation effect not included. 

3. 794 K, 1.0 cm heater with 0.635 cm thick top and side walls, and 1.905 cm thick 
bottom wall. 

4. 1.0 cm heater with 0.635 cm thick top and side walls, and 1.905 cm thick 
bottom wall. Parabolic heater temperature profile (from actual measurements) 
used with 772 I\ centerline temperature (actual radial temperature profile given 
on Figure 7.1). 

Unless otherwise specified above (e.g., simulation #2), all runs included the effects 
of radiation (assuming all heater and wall surface emissivities equal to 1.0), variable 
properties, and wall thermal inertia. A two-dimensional explicit scheme was used 
for the thermal conduction in the walls. The right side of the heater is assumed to 
be insulated from the remainder of the top wall. The outside of the bottom wall is 
assumed to be held at constant temperature (300 K). The outside of the top and side 
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walls exchange heat to the ambient environment via Newton Cooling (this heat loss 
is negligible compared to the heat flux on the inner walls for the relatively short-time 
simulations conducted). The results of these simulations and the experiments which 
they model are presented below. Note that the left side of the computational domain 
corresponds to the axis of symmetry of the container. 

As shown on Figure 7.2, one toroidal cell is initially formed in the container 
for simulation #1. After approximately 4 seconds a second cell begins to form near 
the side wall at near mid-height of the container. This cell eventually splits the 
flow pattern into three separate toroidal-shaped vortices. Figure 7.2 shows a velocity 
vector representation of the flow fields, with all magnitudes of the vectors equal in 
order to emphasize the qualitative pattern of the flow. In actuality the velocities just 
below the heater are on the order of 3 cm/s while the velocities in the middle and 
lower vortices are on the order of 0.2 cm/s and 0.05 cm/s, respectively. Therefore, in 
order to get a better quantitative and qualitative picture of the flow field, it is helpful 
to plot contours of the nondimensional stream function, defined by 

rp = — J prvdr , ip = r J pudy (7.1) 

Plots of this type indicate that the maximum value of the stream function in the 
middle vortex increases slightly until t = 15 seconds, decreases between t = 15 seconds 
and t = 50 seconds, and then increases steadily for times greater than 50 seconds. This 
is shown on Figures 7.3 and 7.4. It appears that for t < 50 seconds, the middle vortex 
is partially entrained by the upper and lower vortices. However for t > 50 seconds, 
the middle vortex strengthens and the amount of entrainment decreases. Notice 
that at t — 10 seconds the ratio of the local maxima of the stream function for 
the three vortices (upper:middle:lower) is 126:29:5 (approximately 25:6:1), while at 
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Figure 7.2: Flow patterns from simulation #1 for (a) t = 0.5. (b) t = 4, (c) t = 15, and 
(d) t = 50 seconds. All velocity vectors shown are made equal length to emphasize 
the qualitative pattern of the flow. 
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t = 50 seconds this ratio is 146:16:48 (approximately 9:1:3), and at t = 200 seconds 
this ratio is 176:39:53. 

As shown on Figures 7.5 and 7.6, only one large toroidal cell is distinguishable 
during the first approximately 50 seconds in experiments using the 682 K heater (note 
that this is the centerline temperature of the heater; the exact radial temperature 
profile of the heater has not been measured). At later times, a second vortex appears 
in the lower half of the container. This vortex is much weaker than the upper cell, 
and in some experiments the upper and lower vortices were formed inside one large, 
“figure-8” -shaped cell. The general trend of the experiments performed to date with 
the 1.0 cm and 1.2 cm heaters and the 10 cm (diameter and height) container has been 
that at relatively low heater temperatures (below approximately 650 K) only one large 
toroidal cell forms in the container. Between 650 K and 750 K the single cell starts 
to form a “figure-8” pattern, and for heater temperatures between approximately 
750 K and 820 K two toroidal cells separated by a saddle point appear inside a 
large “figure-8” -shaped cell. For the most part experimental results have shown good 
reproducibility. However, some experiments indicate that soot formation on the inner 
walls of the container (due to the smoke particles used for flow visualization) may 
appreciably affect the radiation heat transfer. We would expect that as the Grashof 
number increases (e.g., by increasing the heater temperature) or as the radiation 
heat transfer increases for a given Grashof number, the number of recirculation cells 
should increase. Over the temperature range investigated the number of recirculation 
cells found experimentally seemed to vary from one cell to two cells. In most tests 
increasing the heater temperature decreased the height of the upper vortex; however in 
some tests the height of this cell increased (extended much lower toward the bottom of 
the container) compared to runs made at a lower heater temperature. One parameter 
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(c) t — 30 seconds 


(d) t = 50 seconds 


Figure 7.5: Flow patterns obtained in experiments using 682 I\ heater for (a) t = 4, 
(b) t — 15, (c) t = 30, and (d) t = 50 seconds. Only one toroidal-shaped vortex is 
distinguishable. Smoke tracers and low-power laser light sheet were used for flow 
visualization. 
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(a) t = 60 seconds 


(b) t = 80 seconds 


Figure 7.6: Flow patterns obtained in experiments using 682 K heater for (a) t = 60, 
(b) t = 80, (c) t = 100, and (d) t = 200 seconds. Outline of relatively weak lower 
vortex is faintly distinguishable. 


(c) t = 100 seconds 


(d) t — 200 seconds 
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that may significantly affect the flow pattern may be the emissivity of the heater 
and wall surfaces. These values may change over the course of the experiments from 
an “ultrasonicallv cleaned” value to a soot-covered value due to the use of smoke 
particles for flow visualization. 

In order to get a lower bound on the effect of radiation on the flow pattern, the 
above simulation was run without including the effect of radiation (simulation #2). 
The result, as shown on Figure 7.7, is that one large toroidal cell forms in the con- 
tainer, the strength of which changes very little from t = 5 seconds to t = 100 seconds. 
A comparison of the nondimensional stream function contours with and without ra- 
diation for t — 100 seconds appears on Figure 7.8. The maximum value of the stream 
function for the single vortex without radiation is only 8% less than the local max- 
imum value of the upper vortex for the case including radiation. However, the flow 
pattern is qualitatively very different since only only large cell is formed instead of 
three. Thus the qualitative flow pattern is a strong function of the radiative heat 
transfer. 

In the first simulation described above, three vortices formed with the middle 
vortex being the weakest and smallest of the three. However, in the experiments 
conducted with the 682 I\ heater at most only two vortices formed inside a large 
“figure-8” -shaped cell. Although the radial temperature profile of this heater was not 
measured (due to equipment failure), another experiment was run using a slightly 
hotter heater with a known radial temperature distribution. Two computer simu- 
lations were run to investigate the differences in flow pattern between using a con- 
stant heater temperature profile and using the actual radial temperature distribution. 
Figure 7.1 shows the heater temperature profiles used for these two simulations. The 
794 K heater temperature corresponds to the thermocouple reading of the centerline 
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fX . 


(a) t = 5 seconds 



(b) t = 100 seconds 


Figure 7.7: Nondimensional stream function (tp) contours from simulation #2 (radi- 
ation not included) for (a) t = 5 seconds and (b) t = 100 seconds. One large toroidal 
cell forms in the container, the strength of which changes very little throughout the 
simulation. 


77 





temperature of the heater while it was in place during the experiment. An optical 
radiometer was then used to measure the radial temperature profile of the heater 
after the heater and water-cooled jacket were removed from the top wall of the test 
cell. Using this measurement technique, the centerline temperature of the heater was 
found to be 772 K. 

Using the constant-temperature heater profile, the flow pattern developed in a 
qualitatively identical manner as the previously reported simulation with the 682 K 
1.0 cm heater. However, when the radially- varying heater temperature profile was 
used, a “figure-8”-shaped flow pattern resulted with a saddle point separating two 
clockwise rotating vortices. Figure 7.9 shows a comparison of the two flow fields 
after t — 60 seconds in terms of the nondimensional stream function (note that neg- 
ative values of the stream function correspond to clockwise rotation). It was found 
that, qualitatively, the flow patterns did not change significantly after t = 20 seconds. 
Although no photographs were taken of the experiment with this heater temperature, 
hand-drawn sketches of the flow pattern and reports by the experimenter confirm that 
a “figure-8”-shaped flow pattern was indeed formed, and that this configuration was 
stable over the (several minute) duration of the test. Therefore based on the ex- 
perimental results and the results of the two related computer simulations, one may 
conclude that significant differences in the predicted flow pattern result depending on 
the assumption of the heater temperature profile for a given centerline temperature. 
In order to obtain an accurate comparison of the experimental results, the actual 
radial heater temperature distribution should be obtained. 
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Chapter 8 


Conclusions 


The principal objective of the experimental effort is to validate a computa- 
tional model which simulates transient heat transfer and fluid flow in an enclosure 
containing a single-phase or two-layer gas-and-liquid system heated unevenly from 
above. Toward this goal, both minimum and desired requirement levels were defined 
to guide the design and conduct of the experiments. The rationale presented for the 
minimum requirements was typically based on the need for verifying the qualitative 
results (e.g., flow patterns) of the computer code, whereas the desired requirements 
were established in order to verify the quantitative results (e.g., magnitudes of ve- 
locity and temperature) of the code. Specifications were provided for experiments in 
three separate settings: (a) a normal gravity laboratory, (b) the NASA-LeRC drop 
towers, and (c) a space-based laboratory. Most of the experiments will be conducted 
under normal gravity conditions. However, in order to study the quasi-steady state 
behavior of the preignition problem and to validate the computer program over the 
required range of the Grashof number and the inertia of the flow, eight tests have 
been defined for a space-based facility. The 2.2 and 5.18 second NASA-LeRC drop 
towers will be used to obtain transient low-gravity data and to support the space 
flight experimental effort. For example, the required design for maintaining a flat 
liquid surface in microgravity may be developed in the drop tower tests. As discussed 
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in Chapter 3 and Appendix A, most of the parametric studies will be conducted in 
the ground-based tests to minimize the cost of the spaced-based experimental effort. 

Preliminary results of single-phase normal gravity experiments and simulations 
show fairly good agreement. However some evidence exists which suggests that the 
current computational model should be modified to allow for heater and wall surface 
emissivities less than 1.0. Results also indicate that the radial temperature profile 
of the heater, not just the centerline temperature, must be measured in order to 
accurately simulate the experiment. 
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Appendix A 

Description of Preliminary 
Hardware Concepts 


The ground-based and space-based experiments consist of the following hard- 
ware: test cell, fluids (liquid and gas), heater, flow visualization system, temperature 
measurement system, liquid-gas interface control system, and a data acquisition and 
control system. Hardware concepts for these tests are described in this appendix. 

Test Cell . Previous experimental efforts at NASA-LeRC have utilized plexiglass 
test cells of both planar and cylindrical geometries. In these tests a rectangular tank 
surrounding the test cell was partially filled with an index matching fluid in order 
to minimize distortion due to refraction in the liquid phase. The 1-g tests used 
10 cm inner diameter cylindrical cells of height 10 cm with the wall being 1/8- 
1/4” thick. The rectangular test cell geometry has the advantage that it simplifies 
flow visualization and data reduction. However, the cylindrical geometry has the 
important advantage that an axisymmetric flow is easily achieved within a small 
packaging volume. With the planar geometry, if a rectangular heater strip spans the 
width of the top of the test cell, then the width of the container must be relatively 
large to reduce edge effects. 
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The container must be made of an optically transparent material to allow for 
measurement of the velocity field. Testing in 1-g has shown that Freon 113 (the liquid 
selected for ground-based tests) crazes plexiglass. Hence if Freon 113 (also known as 
R113) is to be used for space-based experimentation, a glass or sapphire container 
should be used. Glass has the disadvantage of being fragile while sapphire has the 
disadvantage of being expensive. Since the preignition tests will not involve flames, 
plexiglass, glass, and sapphire all meet the temperature stability requirements. 

Packaging constraints severely limit the size of the test cell for microgravity 
experimentation. The cylindrical capsule of the 5.18 second drop tower limits the ex- 
periment volume to 1 m diameter and 1 .5 m height. The inside dimensions of a Shuttle 
“special can” are approximately 31” (0.79 m) diameter and 57” (1.45 m) height. The 
minimum dimensions of the test cell are governed primarily by flow visualization re- 
quirements. Tests have also indicated that in microgravity experiments the smaller 
the test cell, the easier it is to achieve a quiescent, “flat” interface. Specifically, the 
required reduced gravity time to reach an equilibrium shape goes as d 1 5 where d is 
the diameter of the test cell. The selected container dimensions will be determined 
after further study of the component packaging and after further assessment of the 
velocity and temperature measurement systems. 

No special provisions are required for maintaining insulated or constant wall 
temperature boundary conditions, as the computational code accounts for general 
two-dimensional conduction and thermal inertia in the wall. Computer simulations 
have indicated that the temperature of the liquid at the bottom of the container does 
not change more than 0.01°C during the first hour of an experiment. This result was 
obtained for the case of a 10 cm height of liquid R113 with an 800 K heater, run 
for one hour in a simulated 1-g environment. The outer side walls, however, heat up 
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to almost 340 K when a 1/4” thick plexiglass test cell is used in a 300 K ambient 
environment for the above simulation. Heat loss due to natural convection to the 
environment is therefore appreciable in quasi-steady state 1-g tests. Requirements 
for ambient temperature have already been given in this document. 

Fluids . To promote stable, measurable fluid motion and to simplify analysis, 
the selected fluid should have the following characteristics: large dcr/dT, small viscos- 
ity, high decomposition temperature, purity, insensitivity to contamination, and well 
known liquid and vapor properties. It would also be desirable that the liquid have 
a near 90° contact angle with the test cell material, although this is not a require- 
ment since special provisions will be made to “pin” the interface (see “Liquid- Gas 
Interface Control System” section). Of course, the selected fluid must also satisfy 
the usual array of health and safety requirements. For space- based experimentation, 
nonflammability is critical. In other space experiment work, silicon oil is proposed 
for use [13]. This fluid is unsuitable however for the present effort as it is a complex 
mixture of polymers for which little or no vapor property information is available. 

An alternative is the halogenated hydrocarbons (tradename: Freon), hereafter 
referred to as refrigerants. Properties for these fluids are well documented with one 
attractive candidate, R113, being exceptionally well documented [18] and well mod- 
elled [19]. This fluid also has a very high rating for being nonflammable and nontoxic 
[18]. As mentioned previously, ground tests with R113 indicated that this fluid crazes 
the surface of a plexiglass container over a few drop tower tests. This is not a problem 
for ground-based experimentation, since test cells may be readily replaced between 
tests. However, volume limitations on the Shuttle or Space Station make spare test 
cells an impractical option for space-based experimentation. For these tests either 
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R113 will be used with a glass or sapphire container, or else alternative fluids will be 
selected for use with a plexiglass test cell. 

Air will be used as the gas phase in all the experiments. Although some difficul- 
ties arise for flow visualization with a two-phase system, using air for the gas phase 
is most practical on a scientific basis. 

Filling the test cell with liquid to a specified level is simple in 1-g experiments. 
However, additional hardware is required for the reduced gravity experiments. For 
drop tower tests the ability to fill the test cell to a specified liquid level via remote and 
reproducible means is required. At the present time a gravity-feed system has been 
designed but not fabricated or tested. For space- based experiments, a different system 
is required. One concept for such a system may operate as follows: a bladder tank is 
initially filled with the test liquid. A flow meter and controller monitor, measure, and 
control the flow and flow rate. A fine mesh screen is installed near the bottom of the 
test cell. The flow rate of the liquid is sufficiently slow (on order of TBD cc/sec) so 
that the surface tension of the liquid anchors the liquid to the screen and walls of the 
test cell, thus preventing the formation of liquid droplets which could wander off into 
the vapor portion of the container. Such a system has been designed and tested for a 
surface-tension-driven space-based experiment [20]. The liquid-gas interface control 
system prevents the liquid from migrating up the walls past the desired fill level. 

Heater . The experiments use a constant- temperature heater operated at either 
350 K, 575 K, or 800 K. This temperature range was selected to provide a variation of 
one order of magnintude in figLAT when the liquid is initially at 300 K. Fluctuations 
in the heater temperature should not exceed ±1% over a period of less than 1 second 
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(the data sampling period). A constant-temperature heat source is used by the com- 
putational model instead of a constant-heat-flux source because the former is more 
representative of actual flames. Normal gravity tests at NASA-LeRC have used a 
heater of shielded nichrome wire brazed to a stainless steel plate. Ceramic cement 
insulation reduces heat loss to the environment, and a temperature-controlled water 
jacket is used to both prevent the test cell material from reaching a melting temper- 
ature and to concentrate the heating from above. Cylindrical heaters of diameters 
1 cm, 2.5 cm, and 5 cm have been constructed for use in cylindrical test cells. However 
for the rectangular test cell which may also be used in ground-based experiments, a 
rectangular heater strip is required. Such a rectangular heater would be surrounded 
by a water-cooled jacket for the purposes mentioned above. The dimensions of the 
heater and/or the emissivity of the heater will be varied in the experiments in order 
to provide for a variation of the radiation heat transfer in the experiments. 

Flow Visualization System . Requirements for the flow visualization system dif- 
fer slightly between the three experiment settings: 1-g, drop tower, and space-based. 
Certainly the goal of the system is the same for each case: to track the flow of the 
gas and liquid phases in order to obtain velocity measurements over time throughout 
the test cell. As a minimum requirement visualization of the general flow patterns 
in each phase is needed. However, different constraints are involved for the three 
experiment settings. For the space-based experiment, packaging and automation are 
the key concerns. The available volume of a Shuttle “special can” is highly restrictive, 
and operator (crew) time is restricted. The flow visualization system for a Shuttle 
experiment must survive launch and be reliable, as crewmembers should not have 
to make any adjustments on-orbit. Drop tower experiments are also appreciably 
volume-limited. The drop tower imposes two additional requirements to the design 
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of the flow visualization system: the system must survive approximately a 30-g shock 
at the bottom of the drop, and the system must function effectively after at least a 
3 minute 1-g “hands-off” hanging time prior to the drop. The 1-g experiments are 
the least restrictive of the three. They have basically no volume constraint, require 
no automation, and undergo no shock or change of gravitational acceleration. 

The flow visualization systems selected for the three types of experiments each 
consist of four basic components: tracer particles, a method for injecting and stirring 
the tracer particles and droplets, a light sheet for illuminating one plane of the flow, 
and a camera or cameras for recording the flow visualization. A grid should be used 
on the test cell to help quantify the visualization measurements. The design of the 
flow visualization system components is discussed below. 

1. Tracer Particles 

Feasibility tests were performed to determine if tracer particles or droplets 
could be employed for flow visualization. Of particular difficulty was the tracer 
measurement of the gas phase. The common technique of smoke tracers was 
investigated and discarded when it was observed that the smoke contaminated 
the liquid R113 surface. This contamination effectively eliminated the surface- 
tension-driven convection because the surface-tension variation due to tempera- 
ture variation was counter- balanced by that due to the concentration variation. 
Therefore smoke tracers will only be used in single-phase tests and two-phase 
tests designed to obtain results without surface-tension-driven convection. In 
place of smoke tracers, small R113 droplets were injected into the gas phase 
near the stagnation point of the heater surface. These droplets, being the same 
material as the bulk liquid, did not contaminate the liquid surface. They were 
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also sufficiently small, on the order of 1-5 micron, so that they allowed obser- 
vation of flow patterns. Figures A.l and A.2 show photographs taken of video 
tape frames from a 1-g experiment. The toroidal-shaped circulating cell in the 
top half of the gas phase is driven by buoyancy. The “notch” in this cell was 
seen repeatedly in several tests, although it was unpredicted with the current 
simulation. However, later experimentation with the orientation of the light 
sheet indicated that the power of the laser light caused this “notch” in the flow 
field. 

In anticipation of velocity field quantification, an error analysis was per- 
formed to determine how well the tracer droplets could track the vapor stream- 
lines. The largest sources of error were found to be gravitational settling (for 
1-g tests), on the order of 2 mm/sec for the largest droplets (5 micron), and 
thermocapillary-driven motion very near the heater surface. 

The tracer particles selected for flow visualization in the liquid phase are 
1 micron aluminum oxide (A/ 2 0 3 ) particles with about 0.5 mg/100 ml concen- 
tration level. This technique has been used successfully in the 1-g experiments 
which were conducted [21]. In order to minimize distortion due to refraction, a 
rectangular chamber surrounding the test cell was filled with an index matching 
fluid (liquid phase only). Due to diffraction, information on velocities within 
1 mm of the surface is lost. In future ground-based tests, surface velocities will 
be obtained by sprinkling aluminum oxide particles on the surface with obser- 
vation from below. Measurement of liquid surface velocities for the space-based 
experiment is desirable but not required. 
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streaklines 


gas/liquid interface 


Figure A.l: Streaklines formed by R 1 1 3 droplets injected into gas phase of 1-g two- 
phase experiment (vertical streaks are reflections of light off the test cell). Photo 
displays both liquid and gas phases at t = 7 seconds after heater is started. The 
liquid phase has aluminum oxide particles in suspension. 
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Figure A. 2: Closeup photo of one-half the gas phase at t = 25 seconds in 1-g two- 
phase experiment. Streaklines are formed by R113 droplets (the large vertical streaks 
are reflections of light off the test cell). The left side of the photo is the axis of 
symmetry of the test cell. 
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2. Injection of Tracer Particles 

In the 1-g tests performed to date, aluminum oxide particles are injected 
and mixed into the liquid by hand. The vapor tracer droplets are injected by 
hand via a hypodermic needle in the center of the heater. For drop tower and 
space-based tests, a remote injection and stirring capability is required for the 
tracer particles. 

3. Light Sheet 

The laser light sheet is produced by an argon or helium-neon laser, a cylin- 
drical lens, and a converging lens, as shown on Figure A. 3. The thin sheet of 
light (about 1 mm thick) illuminates a single plane of the flow in order to trace 
the particle velocities in the plane. In the 1-g tests two light sheets oriented in 
different planes will be used to check for axisymmetry of the flow. The effec- 
tiveness of the flow visualization system depends on both the intensity of the 
light source (i.e., the laser power) and the ability of the recording system (i.e., 
the camera) to detect low light levels. The laser selected for drop tower experi- 
mentation must survive the 30-g impact of the drop. A commercially available, 
DC-powered helium-neon laser which meets this requirement has been identi- 
fied for drop tower use. The 5 mW laser, including power supply, is 18” long 
with a 1.5” diameter. Because the experimental volume of the Shuttle “special 
can” is comparable to the drop tower package, this laser is a good candidate 
for the Shuttle experiment. Another good candidate is a diode-pumped yag 
rod laser. This laser represents the current state-of-the-art in small, relatively 
high intensity (about 20 mW) lasers which do not require much power. Final 
selection of the laser for space-based experimentation is TBD. 
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Figure A. 3: Schematic view of the test cell arrangement. Light sheet serves to illumi- 
nate single plane of axisymmetric test cell. Not shown is rectangular chamber around 
test cell which serves to reduce distortion due to refraction. 
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4. Recording System 

Within the volume constraints of the experiment, a tradeoff exists between 
laser intensity and the ability of the recording system to detect low levels of 
light. The selection of a low-power (5 mW) laser for drop tower (and perhaps 
Shuttle) experiments results in the requirement of a motion picture camera with 
high resolution at very low light intensities for these tests. To determine flow 
speeds on the order of 10 cm/sec within ±2 mm/sec accuracy, the camera speed 
should be at least 100 frames/sec with exposure times less than approximately 
50 fisec. This maximum exposure time represents the time required for a tracer 
droplet or particle to move 5 /im at a velocity of 10 cm/sec. As previously 
mentioned, liquid surface velocities may not be measurable in drop tower or 
Shuttle experiments. If this is the case, the maximum velocity in the liquid at 
least 1 mm below the surface may be only of order 1 cm/sec. Hence camera 
exposure times of less than 500 /^sec may be required. 

Temperature Measurement System . Thermocouples will be used to obtain tem- 
perature measurements of the ambient air, the heater surface, the container walls, 
and selected locations in the liquid and gas phases. Although thermocouples are “in- 
trusive,” no effect of their presence on the flow pattern was observed in 1-g feasibility 
tests. Because thermocouples provide information at only a few discrete locations, 
interferometry was explored as a technique for providing whole field quantification in 
the 1-g experiments. Interferometry is not planned for drop tower or space-based tests 
due to concerns about optical alignment and instrumentation survivability. Of the 
optical methods available for temperature measurement, holographic interferometry 
is preferred for the 1-g experiments. This technique provides whole field quantifi- 
cation non-intrusively. Additionally, a rainbow schlieren deflectometer was used in 
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feasibility tests to provide qualitative information on the temperature gradient field 
in each phase. In these tests a heated wire was placed in air just above a liquid silicon 
oil surface. A holographic interferometer with electronic heterodyne data reduction 
was also used in these tests to quantify the temperature fields near the wire with 
a spatial resolution of 100 ^m. The measurement accuracies in the gas and liquid 
phases were 0.6% and 1.8%, respectively. The somewhat larger variance in the liquid 
phase resulted from the uncertainty in the determination of refractive index based 
upon temperature, whereas accurate tabulated values were available for the gas phase 
measurements. 

For wall temperature measurements, chromel-constantan, type E, 30 gauge wire 
thermocouples will be used with wires running up the side walls. For liquid and gas 
phase measurements, the wire leads should go in a horizontal plane across the con- 
tainer and through the opposite walls of the test cell. Threaded plugs made from the 
container material (e.g., plexiglass) may be used to install and remove thermocouples 
used for the liquid and gas phases. Fine thermocouple wires (approximately 75 fim 
diameter) with butt-end junctions may be used for the ground-based tests. However, 
this set-up may be inadequate for drop tower and space-based tests. 

The measurement of surface temperature distribution with, say, a thin thermo- 
couple probe would be highly desirable and achievable in 1-g. However, positioning 
of the probe close to the surface is very difficult to do in reduced gravity because the 
exact location of the free surface cannot be specified a priori. Even with the proposed 
liquid-gas interface control system, the free surface may deform 1-2 mm. Moreover, 
the free surface deforms a certain extent due to fluid motion and g-jitter. The probe 


i 
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could be positioned manually after the flow is started, but doing this at several ra- 
dial locations for the surface temperature distribution measurement may be too time 
consuming. For the drop tower experiment, such a manual positioning is impossible. 

Due to the aforementioned problems associated with fluid temperature mea- 
surement, the reduced gravity experiments may only incorporate provisions for wall 
temperature measurements. 

Liquid-Gas Interface Control System . Reduced gravity experiments introduce 
the problem of liquid-gas interface surface curvature. In order to prevent such curva- 
ture and to make the comparison of earth and reduced gravity experiments straight- 
forward, special precautions are required to “pin’' the interface and maintain a flat 
shape. Drop tower experiments have been conducted in the NASA-LeRC 5.18 second 
facility which use internal flanges in the test cell (see Figure A. 4). These experiments 
were successful in significantly reducing surface curvature, as shown on Figures A. 5 
and A. 6. Some surface motion was observed throughout the first four seconds of the 
drop. It is believed that this motion will be eliminated by more exact filling tech- 
niques. The internal flanges will be used in both the drop tower and space-based 
experiments to maintain a flat liquid-gas interface with a 0.04 radian maximum local 
angle of the liquid surface. The positioning of these flanges along the container walls 
should be adjustable to allow for a variation in the liquid fill level. 

Software for Data Collection . Temperature and flow visualization measurement 
data from the experiments must be collected and analyzed. A commercially-available 
data acquisition and control system has been identified for the drop tower experi- 
ments. This system is capable of sampling at the rate of 100 readings/second/channel 
on each of 8 channels. With its digital I/O, it is capable of control functions. The 
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Heater 



Figure A.4: Internal flanges used to “pin” the gas-liquid interface. Peak-to-peak 
length and distance from the wall are 1/16 inch (1.6 mm). 


97 



ORIGINAL pagE JS 




Figure A. 5: (a) and (b) display the liquid position in a 10 cm tank at various times 
in a 5.18 second drop, (a): about 0.1 seconds after drop start, (b): about 4 seconds 
later; note the liquid has crept along sidewall, forming a nearly hemispherical shape. 
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Figure A. 6: (a) and (b) are analogous to (a) and (b) in Figure A. 5 except that internal 
flanges have been installed in the tank. These prevent the liquid migration and serve 
to “pin” the liquid-gas interface nearly flat. 
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selection of the data acquisition and control system for the Shuttle experiments has 
not yet been finalized. 

In order to quantify the flow visualization, software has been developed for use 
with an NEC motion analyzer which digitizes 16 mm motion picture frames. Data 
collected on the Shuttle will be evaluated on the ground. 

Safety Precautions 


Since the experiments will not involve flames or flammable or toxic substances, 
no special handling requirements or safety precautions are required (non-decomposed 
R113 is classified as nonflammable and nontoxic by [18]). Standard OSH A (for 
ground-based) and Shuttle or Space Station (for space- based) guidelines should be 
followed for the operation of the experiments. A temperature control and automatic 
shutoff system should be used in the heater design for overtemperature protection. 
Although the heater will be insulated in the experiment, computer simulations have 
indicated that the outside wall temperature around the gas phase may reach ap- 
proximately 340 K in a one-hour experiment. Therefore care should be taken when 
handling the test cell during and/or after an experiment. 
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Appendix B 

Derivation of Energy Equation 


In order to write the thermal energy equation in terms of excess enthalpy, we 
define for the species (all quantities dimensional unless otherwise specified) 

hi, excess — Cp,(T — T°) = c Pi 7\ (B-l) 

where 

h, = = h-h° 

t 

hi = h° + ft,-, ex cess = h° + h it i 
h° = hi(T°) = formation enthalpy 
T° = T m = 300 K (reference temperature) 
h = ZYihi 

i 

h° = EW 

t 

c, = (yrp) 

c r = (B.2) 

* 

The thermal energy equation will be derived in rectangular coordinates for the 
two-dimensional planar geometry. The results will then be generalized to apply also 
to the cylindrical geometry. 
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To obtain the thermal energy equation, we derive the total energy and mechan- 
ical energy equations, and subtract the latter from the former to get the thermal 
energy equation. The following notation will be used: 


ii t = internal energy per unit mass of species i 


_j_ 2 

9 1 ’* 


r TT — 


T yy ~ 


T y x — 


6 

Qv ~ dV 


= kinetic energy per unit mass of species i 


du 2 - _ 

P' v 


d 
d 

T xy = /' 


dv 2 , 
dy 3 V ‘ t; 
dv du 
dx dy 


P — T x 


— r. 


XV 


— T 


yx 


P - T 


yy 


volumetric heat source term 

(energy production per unit volume per unit time) 


$ = viscous dissipation function (see [16], eqn A. 16) 


Notice that the sign convention for the shear stresses is opposite to [15]. For a two- 
dimensional control volume dx ■ dy , the conservation of total energy is 


> 

rate of 


' > 

net rate of 



accumulation 


internal and 


net rate of 

< of internal 

► = 4 

kinetic energy 

> + < 

heat addition 

and kinetic 


in by convection 


by conduction 

energy 


and diffusion 

< > 


t j 
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+ < 





net rate of 


net rate of 

work done on 

► + < 

heat addition 

element by 


by sources 

surroundings 

t * 


^ d 


(B.3) 


An Eulerian point of view will be used for the derivation (fluid flows through a control 
volume which is fixed in space). 


The rate of accumulation of internal and kinetic energy is 


d_ 

dt 


+ \ v l) Y ip 


dx dy 


(B.4) 


where p = /Vnjv = density of the mixture and (following [22]) 

= Vi’Vi 

Vi = v + V t 
V = Yl Y iVi = VA 


Vi = 


Yi 


(B.5) 

(B.6) 

(B.7) 

(B.8) 


In equation (B.6), v, is the local average velocity of species t (relative to stationary 
coordinate axes), v' is the mass-weighted average velocity of the fluid mixture [given 
by equation (B.7)], and K is the diffusion velocity of species i [given by equation (B.8) 
which follows from Fick’s law of diffusion]. 

Substituting equations (B.5)-(B.8) into equation (B.4) gives 


d_ 

dt 


Y^{iii + )f)Yip dxdy = dx d V + 


d_ 

dt 


ip(vV)£y; + P YiW-Vi)Yi + 


dx dy 
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For diffusion in a nonstationarv medium, the total mass flux of species i (denoted 
by n,) is given by 

n, = pY,v' - pV^ = pY(v' + Vi) (B.13) 

The first term on the right hand side of equation (B.13) is the convection of the mass- 
weighted average fluid flow, and the second term is the diffusive contribution which 
represents the mass flux of the i in species relative to the mass average velocity v . The 
amount of internal and kinetic energy which flows into the Eulerian control volume 
through the z-face by convection and diffusion is 

d y 

l J x 

while the amount of energy which flows out through the x + dx face by convection 

and diffusion is 

+ + £»?>] ^ 

* -• x+dx 

Hence for the ar-direction, the net rate of energy in by convection and diffusion is 

Y.P Y ^ v 'x + K>)(«. + ^v, 2 ) dx dy 
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Similarly, the net rate of energy in by convection and diffusion in the ^/-direction is 


d_ 

dy 


2Z^K + k.»)(«, + ^, 2 ) 


dx dy 


Therefore the total net rate of energy in by convection and diffusion is 


-V- 


+ ?<)(* + y*?) 


dx dy 


In the above equation vf = where u,- is given by equation (B.6). Therefore, 
substituting equation (B.6) into the above equation yields 


-y. 




dx dy = 


= -y. 

= -v« 


pE Y iW + ^)(«. + + v'-Vi + ^Vi-Vi) 


dx dy 


pv' ]T Yi^i + p E YiUiVi 


dx dy 


_y. 


1 


52 if + pv 52( v ''Vi) Y i 


+ 5 pv'ZW-vm 


dx dy 


_y. 


\ p(v'-v') 52 y.Vi + p 52 yy&'-Vi) 




dx dy 


Using equations (B.5), (B.6), and (B.9)-(B.ll), the total net rate of energy in by 
convection and diffusion is 


-V- (pv‘ u + pE dx dy 


-V* 

_y. 


i p v\v > -v > )^\pv'E,{VrV i )Y i 


dx dy 


(B.14) 


p 52 yMv'-v ; ) + ip 52 y<Vi(v,-vd 


dx dy 
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Denoting 


k = the mixture thermal conductivity, 
the rate of heat addition in by conduction is 

(-*f) **(-•*)/■ 

while the rate of heat addition out by conduction is 

(-*£) *+(-*£) <** 

V 9xJ x+dx V d yJy + dy 

Therefore the net rate of heat addition by conduction is 


(V-jtVT)dx dy 


(B.15) 


Notice that radiative heating of the gas has not been included in the analysis. 


Work is done on the element by body forces (gravity) and surface forces (pressure 
and shear forces) 1 . The rate at which work is done by the body forces is the force 
times the velocity: 

pg'v'dx dy 

Similarly, the rate at which work is done on the element at a face is the vector dot 
product of velocity and force. For example, on the x-face of the element the force 
in the x-direction is (P — t xx ) x i. Therefore the net work done on the element at the 
x-face in the x-direction is 

dy{P - T xx ) x i m v' = v' x (P - T xx ) x dy 

Similarly, work done on the element at the x-face in the y-direction is 

{-T xy dy) x j-v' = -v' y r xy \ x dy 

l \n the derivation of the work done on the element, the contribution of work done by diffusion 
has been neglected to simplify the analysis. 
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Therefore, for the x-face the work done on the element is 

( v'*F x ) x dy 

The convention for positive shear stresses is such that a positive stress has 
opposite direction in opposing planes. This is usually the case. However in some 
instances, such as if the center of the control volume is along the centerline of a 
Poiseuille flow, the direction of positive shear is the same for opposite planes of the 
control volume. To be precise in general, therefore, the net force from pressure and 
shear on the x and x + dx faces is of order dy. Using the previously defined positive 
shear directions, the work done on the control volume at the x + dx face is 

~(v'‘F x ) x+ d x dy 

The net work done on the element at the x and x + dx faces is therefore 

-~(v'-F x )dxdy 

For the y and y + dy faces, the total work done on the element is 

d - - 

-—(v’'F y )dxdy 

Oy 

By adding the contributions at all four faces, we get the total work done on the 
element by surface forces: 

- + -^{v'Fy) dxdy 

Therefore the net rate of work done on the element by the surroundings (body forces 
and surfaces forces) is 

| P9'v - ^(v'-F x ) + -^(v'‘F y ) ^jdxdy (B.16) 
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The net rate of heat addition by sources is given by the volume of the element 
multiplied by the volumetric heat source term, Qv- 


Q v dx dy 


(B.17) 


By substituting equations (B.12)-(B.17) into equation (B.3) and dividing both 
sides by dx dy , we are left with the total energy equation (before applying continuity): 

f tipi + l^ ) + ll^y,Vrv) = 

= —V* (puv + p^2Y,UiVi\ — V* 


^pv'(v'-v') 


+ \pv'Y. y W-V*) 


_V 


P^YYiv'-V,) + -pY:y t Vi(VrV,) 


-fV'kVT + pg'v' — 


Q — - ()-.-> 


+ Q\ 


In order to simplify the equation, we will use the overall continuity equation, 

(B.18) 


+ v'-Vp + pV-S = ^ + V-(pv') = 0 
at at 


The first term on the left hand side and the first term on the right hand side of the 
total energy equation can be combined with the help of equation (B.18) to- give 


J t (pu) + V-(puv') = 

= ul^ + V-(pv‘) 
du „ 

= „l_ + „. vu 


du 

+ p-z 7 + pv *Vu 
c it 


= P 


Du 

Dt 
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Similarly, we can write 


^-(ipv'.v') + V-[^pv'{v'-v')] = 




| + V.(^) 


1 

+ 2 P 


d(v-v) 

dt 


+ u'*V(u'*i/) 


D ,1 -5 - 5 . 

= P*7T~( -v *v ) 

H DV 2 ’ 


and 






| + v. ( ^) 


1 

+ 5 ' 


-EW-K + .'-V^W-V, 

. l x i 

44 

The above equations were each simplified using the overall continuity equation [equa- 
tion (B.18)]. Substituting these three equations, the total energy equation becomes 

D . „ l - *, 1 X?/ r-s -* T -t\ 

/, Di {u + ? ' v) + 2 "Di (EW-KJ = 

= -v- |/.E[>i“.K + YiViW-Vi) + jKiK(K-V;)]| 

+V-kVT + pg-v’ - ^W-F*) + ^('"'-F v )] + Qv (B.19) 


The mechanical energy equation is the dot product of velocity (neglecting me- 
chanical work by diffusion, we use the mass-weighted average velocity of the fluid 
mixture) with the momentum equation [equation (6.5)]. The momentum equation 

(where we now replace V with v‘) is 

— * 

*4 = ~ vp + 1^1 + « (B-20) 
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Substituting the stress tensor into the momentum equation, equation (B.20) becomes 


,£i = _vp + /aT ” 5r ' 


Dt ' \ dx dy , 

In terms of F x and F y , the above equation is 

Dv' df x dF r 


" + ^|7+l^ + ^7 +« 


dx dy 


Dt 


dx dy 


+ P9 


Therefore the mechanical energy equation is 


1 Dv 

pv — = — u 




+ pv<7 


(B.21) 


To get the desired thermal energy equation, we subtract the mechanical energy 
equation [(B.21)] from the total energy equation [(B. 19)] - To do this, it is easiest to 
re-write the contribution of net work done on the element in the total energy equation 
first. Applying the chain rule to this part of equation (B.21), we get 


= — V 




, _ 0 dv' ~ dv 

+ pv ~ F -' & " F »“ a. y 


Using the above identity, subtracting equation (B.21) from equation (B.19) leaves the 
thermal energy equation 

D , „ 1 Dv 1 D (r-y - 

p m (u + 2 V ' “ pv "m + 2 p m V? ' ’ V = 

= -V- L Yimv, + WW-Vi) + iy,K(K- K)]J 

+V-I.-VT - F x ~ - Fy~ + Qv 

We can also substitute for the viscous dissipation, /i$, which is (for two dimensions) 


- - - dv' -* dv' 

lli -PVv=-F,~-F,~ 


(B.22) 
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and also note 


1 Dv_ D 1 -* - 

• dT = , ’5I ( 2 t ’ > 


so that the thermal energy equation in terms of the mixture internal energy becomes 

Du l D ~ A 

p ~m + 2 p m = 

= -v- jpDwtf- + YiViiv'-Vi) + ^(tf-tf)]} 

+V‘fcVr-PV*t? + ^$ + Q v (B.23) 

We can re-write the thermal energy equation in terms of the specific enthalpy by 
using the relationships 

Yihi = YiUi + j (B.24) 

and 


h = '£Y i h i = ii + - 
i P 

Substituting equations (B.24) and (B.25) into equation (B.23) gives 


(B.25) 


P Dt P Dt 


= -V* 


V p , 


- v- + |v;k(k.k)]J 


+V-kVT — PV‘v + p$ + Qv 


Noting that 


Dim = _DP+FD l= _DPP_la t , \ 

Dt Dt p Dt Dt p \dt P ) 


and from the overall continuity equation [equation (B.18)] 

dp , X-, & 1 

— + u *V/> = -pV'v 


we get 


D(P/p) DP - -5 
~P — — — = — — - Pv-v 


Dt 


Dt 
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so that the thermal energy equation simplifies to 

-v- - v- (t?-3) + |v;k(k*k)]| 


+^ + v-wr + ^ + Qv 


(B.26) 


Notice that in the absence of diffusion (V^ = 0 for all i) the thermal energy equation 
reduces to 

D u _» -* 

P — + PV-v' = V-kVT + /iO + Q v 

or 

4=f + ^ vr+ " 4+ ^ 

as we would expect. These two equations apply to a single-component system. 


In order to write equation (B.26) in terms of the excess enthalpy, we substitute 
equations (B.l) and (B.2) to get 




-v- 


p^2v.y,(h° + h u ) 

i 

f + v-fcvr + //$ + q v 


- v- {pTjyy.w-v .) + jwtv’.-K)]} 


This can be expressed in a slightly different form using the continuity equation by 
writing 


= Jt {ph ° ] + - h ° 


| + v. ( „7) 


Since the last term in the above equation is identically zero by the overall continuity 
equation, we get 
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Similarly, 


Dh\ d , , x a * \ 

Substituting the above two equations into our thermal energy equation gives (after 
re-arranging terms) 

J t (ph,) + V-(-kVT + pv'hi) = V- £ PiV, 

- + v-(p?h° + pY, VMK ) 

-V- (,£ VM.hu) + ^ (e W-k) 

-V- ^ + Qv 

Consider the first three terms in brackets on the right hand side of the above equation. 
Since 

h° = constant for each i 
and h° = Y{h° 

we can write 

|(^»)+v- Uh'+p-Ewhi) = 

= E {*? [!</>«) + 1 

= E{*?[|(^) + V.n-]} 

where rT, is the total mass flux given previously in equation (B.13). From the conser- 
vation of species [equation (6.28)], 

r i = oi^pYi) + v * n »' 
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where 


r, = rate of production of the i 
component [kg/(m 3 -s)| 


:th 


we end up with 


-( ph °) + V-\pv'h° + p 


Euw?) =- 


E *fr" 


If we define the dimensional quantity 


^ h°r = heat released by chemical reaction 
* (calculated at temperature T°) 

then the thermal energy equation becomes (after some rearrangement) 

j t (phi) + V-(-kVT + pv'hi) = ^ + <?ch 

V- (/. E feu + p* + Qv - \p£- t (E W-k) 
-v-|/>E[ y .K("'-^) + 5^(^- 17 .)]} + v-£p,k 

In order to use = k' /c' p and <j> = c p T[ = h\ = ^excess* we can write 


V-(-kVT) = V-(-kVTt) = -V- 
= V- 


fc fcr, __ 

— V(cp7\) - ^-Vcp 


L C P 


k„, 1 


r ArT’i 1 

V/i, 

4- V- 

T^VCp 

c p J 


C P J 


and use 


Tl = — = 


(B.27) 


(B.28) 


(B.29) 


P t; pD i; 

so that by substituting equations (B.l), (B.8) and the above two equations into 
equation (B.28), we get 
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+v- Lx>„r, JfoVV',) - LeVZ r | 

-V- j/> £[*#(.?•{$) + i KiKfK-K)]} + v- £ P.K (B.30) 

All variables in equation (B.30) are dimensioned. Note that we have assumed the 
binary diffusion coefficients, T>ij, are equal for all species. If we denote 

£ = %■ + <4 + {^; [d 4 w-) - ie Vc;] } 

then equation (B.30) becomes 

^r(/> *,) + V. (-=rV A, + p = 5* 

The mass- weighted average velocity of the fluid mixture [ 1 / (non-dimensional) or 
(dimensional)] is the same velocity as is used in the other conservation equations. 
Thus the above equation written in standard form [see equation (6.1)] using the 
mass-weighted average velocity components is 

W {p Al) + 7a? [ r (/ M ‘*'-p7j] + 


d ( , k'dh'A c , 

+ w{ pvl "-%w =s * 

To make equation (B.32) dimensionless, multiply it by 

t * 

P»Cp t Tm p*UmCp,T+ 


(B.32) 
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and use the nondimensional variables 


/>! = 



h. c p ,T . 


and 


The result is 



d_ 

dt 


, , , Id 
{/lhl) + 7Tx 



pjf. IzJCp dhA 
U.L. n . dx) 




P_J_p. k'fc, dhA 
U.L. p. dy ) 


tl s' 

P*c p .T.* 


Recognizing that 


pjp- _ 1 

U.L . “ Re. 


the above equation is in the general form of equation (6.3), 

r* d<t>' 


d{p<i>) 1 d 

dt r dx 


d 


H' 


pu4> — 

r* d<f>\ 


+ dy Re. dy 


Re. dx ) 
t. 

P* c p.T.‘ 


with 


<f> = hi — 


fh 

Cp.T. 


and 


r* = 


k'/c p 


(B.33) 


(B.34) 
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Appendix C 


Derivation of Source Term For 
Species Equation With One- Step 
Chemical Reaction 


In the current version of the computational model, chemical reactions are not 
included, and therefore the source term given in equation (6.31) is identically zero. 
That is, 


S* = 



= 0 


However, in future versions of the code incorporating chemical reactions, the source 
term for each species may be developed as follows. Consider a simple one-step chem- 
ical reaction 

npF + uqO npP (C.l) 


where 


F,0,P 

np , no, Tip 

k' 


n = np + no 


= fuel, oxidant, products species 
= moles of fuel, oxidant, and products in reaction 
= specific rate constant 
= order of reaction 
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Then the rate of production of fuel [units: kg/(m 3 s)j is 

(’>)' - * 


d ^hl' F = -n F k'[F] n r[0] no Al' F 


(C.2) 


where M f is the dimensional molecular weight of the fuel and the single brackets 
(e.g., [F]) denotes concentration, C' % . The concentration C\ is given by 


C [ = Pm^i = (P.Y.\ (PY 


M'i 


M. ) \ Mi , 


(C.3) 


where Y. = 1. 


The specific rate constant is given by the Arrhenius expression 

k' = A' exp(-E'jn'T') (C.4) 

where A' is the pre-exponential factor (representing the number of collisons per sec- 
ond), E' a is the activation energy, 7Z' is the universal gas constant, and T' is the 
absolute temperature. 

Substituting equations (0.3) and (C.4) into equation (C.2) gives 

(->)' = -n F p:Ml-M F X exp(-E'jK'T') ° (C.5) 

Similarly, for the oxidant we get 

(r'o)' = -nofZMl-'MoA' exp(-E'Jfc'T') (^) " (^) ° (C.6) 

The reaction rate for the products is 

(>>)' = n pp"Ml- n MpA' exp(—E a /7Z'T ) ' {S&j ° (C.T) 

As can be seen from equations (C.5)-(C.7), the reaction rates for fuel, oxidant, and 
products differ only by a pre-exponential factor. 
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In the nondimensional source term from the energy equation, we had [equa- 
tions (6.48) and (B.27)]: 

p.c r .T. P.Cp.T. i P- ) 


From this one recognizes that we should write 


h° ... (>>) <■ 

K = — = — 


Now if we define 


(C.8) 


PREXP = A' p"~ l Kll~ n t m and EART = -E'JlZT. 


(C.9) 


then the species source term equations become 


t. / nn ^ vn /EART\ ( 

S^f = — (r F ) = -n F Mp- PREXP* exp -jr ■ j ^ 


/EART 

St# = -n 0 M 0 - PREXP- exp (— 

« /EART\ (pY F \ 

S<ti,p = n P Mp*PREXP* exp - — J ^-r-J 


pYA nF ( pYo \ 
Mp) \M 0 ) 

\ (p}V\ nF (pyoY° 

)\Mp) \Mo) 


no 


nr 


( P Yo\ 

\Mo) 


no 


These can be put into a common form 

r, / EART \ / P Yp\ " F ( P Y 0 \ 

S ‘" - S.-PREXP. exp (— ) ( jjr) (^-) 

where 

A A, * 

1 3f = —npMf, Bq = — noAfo, and Bp = npMp 


(C.10) 


(C.ll) 


These coefficients f?, may be thought of as the mass of species i produced in the 
reaction described by equation (C.l). As described previously, these source term 
equations may be used in the Q c h contribution to the energy equation source term 
when chemical reactions are included in the model. 
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